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FOREWORD 


This  counter  pro^rnn  is  ons  of  a  ssriss  of  dipit?.!  cc"”«tcr  "roprisns 
developed  as  in-house  effort  in  support  of  tne  Space  Vehicle  Thermal  and 
Atmospheric  Control  Study*  The  study  is  sponsored  by  the  Flight  Dynamics 
Laboratory  of  the  Research  and  Technology  Division  under  Contract  AF33(657)- 
S9S5  and  is  under  the  direction  of  W.  Uhl  of  the  Enviror.mertal  Control  Branch* 
R*  F..  Sexton  of  SSJD  served  as  Project  Manager  of  the  studv  program.  H,  b. 


r.:  1 


FDL  TOR. 


The  program  described  in  this  renort  represents  the  second  stage  in  the 
development  of  a  general  configuration  factor  computer  program.  This  report 
partly  incorporates  SID  62-393  (ASD  TN  61-101),  which  describes  CONFAC  If 
the  first  program  developed  under  the  Space  Vehicle  Thermal  and  Atmospheric 
Control  Study. 

This  renort  may  also  be  identified  b/  Contractor’s  Report  No,  SID  63-1397, 


ABSTRACT 


A  simple  numerical  method  is  derived  for  the  determination  of  the  geometric 
radiant-interchange  factors  used  in  radiant  heat  transfer  and  illumination,  A 
FORTRAN  II  digital  computer  program  utilizing  this  method  is  developed  which 
provides  a  rapid  and  accurate  moans  of  computation  of  configuration  and  form 
factors.  The  source  of  flux  may  be  any  general  plane  polygon  and  the  receiver 
may  be  any  general  olane  or  nonplanar  polygon,  the  surface  of  an  arbitrary 
polyhedron,  or  an  arbitrary  combination  of  such  surfaces. 

It  is  therefore  possible  to  accurately  determine  configuration  and  form 
factors  from  a  plane  surface  to  another  surface  occluded  by  complex  interven¬ 
ing  surfaces.  Form  factors  are  computed  rapidly  —  averaging  less  than  two 
seconds  on  the  IBM  7094  for  simple,  unobstructed  plane  surfaces,  and  less  than 
30  seconds  for  simple  polyhedra.  Also,  means  are  provided  to  internally 
generate  a  variety  of  regular  polygons  or  poiynedia  and  to  transform  surface 
spatial  coordinates  for  convenience  of  data  entry  and/or  motion  simulation. 
Simplicity  of  data  entry,  flexibility  of  application,  and  economy  of  operation 
are  principal  features  of  this  program.  Sample  problems  illustrating  these 
important  aspects  are  provided. 

This  report  has  been  reviewed  and  is  appi'oved. 


~  BAKER 
"Asst,  for  R&T 
Vehicle  Equipment  division 
AT  Flight  Dynamics  Laboratory 
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d^  Differential  area 

dA-A  From  a  differential  area  to  an  area 
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SECTION  I 


INTRODUCTION 


The  geometric  form  factor,  *s  defined  as  t^e  fraction  of  radiant 

energy  emanating  from  finite  surface  which  is  intercepted  by  another 
surface  A2 . 


$  Flux  received  by  finite  surface  A2 
12  "  Flux  emitted  by  finite  surface  Aj  ^ 

The  geometric  configuration  factor,  ~yz*  defined  in  a  similar  manner, 
except  that  the  emitting  surface  is  infinitesimal,  (sometimes  referred  to 
as  the  point  configuration  factor). 

Flux  received  by  finite  surface  A~ 

C  j,  .  -  -  .  ■  . .  - - — — £a  .... —  <2 > 

Flux  emitted  by  infinitesimal  surface  dA]  v  * 

The  subscripts  denote  the  direction  of  flow  of  net  flux;  c12  and  fj 2  pertain 
respectively  to  the  configuration  and  form  factor  from  surface  A|  to  surface 
A2«  It  is  assumed  that  each  surface  is  isothermal  and  radiates  diffusely, 
i,c.,  follows  Lambert’s  cosine  distribution  law. 

The  "closed-form”  determination  of  the  configuration  or  form  factor  by 
classical  integration  techniques  is  impossible  or  impractical  in  most  situa¬ 
tions,  Experimental  techniques  and  devices  have  been  reported  in  the  litera¬ 
ture  (Reference  1),  and  probably  the  most  useful  is  Pleijel’s  Globoscopc 
(Reference  4),  Experimental  techniques  produce  only  the  configuration  facto! 
however.  Nonetheless,  they  are  useful  for  many  applications  where  only  one 
or  just  a  few  configuration  factors  are  required  and  nominal  accuracy  is 
sufficient. 

However,  if  a  large  number  of  form  factors  are  required  in  a  short 
period  of  tine,  experimental  techniques  are  not  practical.  This  report  pre¬ 
sents  a  numerical  method  and  a  computer  program  which  enables  rapid  and 
accurate  computation  of  configuration  and  form  factors  between  plane  sur¬ 
faces,  and  plane  or  solid  surfaces.  The  source  (surface  1)  may  be  any  general 
plane  polygon;  the  receiver  (surface  2)  may  be  any  arbitrarily  oriented 
general  plane  or  nonplanar  polygon,  the  surface  of  an  arbitrary  solid,  or  an 
arbitrary  combination  of  planes,  nonplanes,  or  solids.  Form  facurs  (which 
nominally  are  derived  from  025  configuration  factors)  are  commuted  rapidly, 
averaging  less  than  2  seconds  by  IBM  7094  rime  ^or  simple  plane  surfaces, 
and  less  than  30  seconds  from  simple  plane  sur> ^ces  to  simple  solids.  Table  1 
compares  solutions  obtained  by  CON FAC  II  to  ti.vse  given  in  Reference  1, 

Manuscript  released  bv  the  author  January  1964  for  publication  as  an  FDL  tech¬ 
nical  Documentary  Report 
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Table  1,  Comparison  of  Configuration  and  Form  Factors  Computed 
by  CONFACII  to  Those  Given  in  Reference  1 


Configuration 
P-1,  X  =  1,  V  =  1 

x  =  o.i,  y  =  o.i 

x  =  i,  y  =  * 

X  =  0.1,  Y  =  0.4 

X  =  1,  Y  =  <®  * 

X  r  o.i,  Y  *  »  * 


P-2,  0 

0 


0 

0 

0 

0 

0 


30°,  L  =  0,  N  =  1 
30°,  l  -  1,  N  =  1 
30°,  L  =  0,  N  =  4 
30°,  L  =  4,  N  =  4 
120°,  L  =  0,  N  =  1 
120°,  L  =  1,  N  --  1 
120°,  L  s  0,  w  s  4 
120°,  L  =  4,  N  =  4 


**P-8,  D  =  4,  L  =  2 

D  =  2,  L  =  4 


X  =  0.1,  Y  =  0.1 

X  =  1,  Y  =  4 

X  =  0.1,  Y  =  0.4 

X  =  1,  Y  -  «  * 

X  =  0.1,  Y  =  »  * 

A-2,  0  =  30°,  L  *  1,  N  a  1 

0  =  30%  L  a  4,  N  «  4 

0  =  120%  L  s  1,  N  =  1 

0  »  120%  L  =  4,  N  =  4 


Computer  (Trapezodial  Rule) 
Reference  1  24  x  24  grid  60  x  60  grid 


0.13853 

0.00314 

0.17527 

0,01147 

0,17678 

0.02488 

0.4665 

0.1759 

0.4665 

0.0964 

0.125 

0.0236 

0.125 

0.0077 

0.276 

0.438 

0.500 

0.724 

0.800 

0.08074 

0.24774 

0.19982 

0,00316 

0.34596 

0.01207 

0.41421 

0.04988 

0.6202+ 

0.3961+ 

0.0870+ 

0.0433+ 


0.138532 

0.003141 

0.17527C 

0.011471 

0.176777 

0.024876 

0.466506 

0.175923 

0.466506 

0.096447 

0.125000 

0.023554 

0.125000 

0.007683 

0.275 

0.436 

0.498 

0.722 

0.799 

0.08055 

0.2472 

0,19972 
0,00316 
0,34559 
0, 01207 
0.40549 
0.04884 

0.61769 

0.39431 

0.08665 

0,04272 


0.138532 

0.003141 

0.175270 

0.011471 

0.176777 

0.C24376 

0.466506 

0.175923 

0.466506 

0.096447 

0.125000 

0.023554 

0.125000 

0,00768* 


0.19981 

0.00316 

0.34590 

0.01207 

0.41075 

0.04946 

0,61878 

0,39450 

0,08662 

0.04235 


S 

*  10  was  assumed  to  approximate  *»  for  computer  run 

**  32  sided  regular  polygon  used  to  simulate  circular  cross-section 

+  These  values  were  obtained  by  numerical  integration  across  surface  Aj, 
according  to  Reference  1 
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The  FORTRAN  II  Computer  Program  described  herein ,  CONFAC  II,  is  a 
follow-on  development  of  an  earlier  version,  C0NF4C  I  (Refer^^e  5).  The 
original  program  has  been  extensively  modified  and  significant  improvements 
in  the  flexibility  of  application  have  been  achieved-  CONFAC  I  was  developed 
principally  to  compute  geometric  form  factors  between  plane  surfaces;  no 
application  to  nonplanar  surfaces  or  bodies  was  originally  intended*  However, 
because  of  the  particular  analytical  apnroach  utilized  and  the  data  handling 
techniques  developed,  it  was  possible  to  use  the  basic  plane-to-plane  program 
to  compute  factors  to  nonplanar  surfaces,  provided  proper  restrictions  were 
observed* 


The  principal  similarities  and  differences  between  CONFAC  I  and  CONFAC  II 
are  as  follows: 

1.  Both  CONFAC  I  and  CONFAC  II  require  that  Surface  1  be  a  plane 
polygon;  it  may  be  arbitrarily  oriented  in  the  coordinate  system 
in  which  it  is  described  (entered)  in  data. 

2*  Both  CONFAC  I  and  CONFAC  II  specify  that  if  Surface  2  is  a  plane 
polygon,  it  may  be  arbitrarily  oriented  with  xcspect  to  Surface  1 
and  within  its  own  coordinate  system. 

3.  Both  CONFAC  I  and  CONFAC  II  require  that,  if  Surface  2  is  a  non¬ 
planar  surface,  then  the  surface  boundaries  must  present  a  valid 
silhouette  from  any  point  on  the  active  side  of  Surface  1. 

4.  CONFAC  I  specifier,  that  no  part  of  nonplanar  Surface  2  may  lie 
below  the  "horizon"  of  Surface  1  when  viewed  from  the  active  side 
of  Surface  1.  CONFAC  II  does  not  require  that  all  of  a  nonplanar 
Surface  2  appeal  above  the  horizon  of  Surface  1,  CONFAC  II  will 
automatically  bisect  a  nonplanar  Surface  2  and  compute  the  factor 
to  only  the  part  which  Surface  1  "sees." 

5.  CONFAC  I  cannot,  in  general,  be  used  to  compute  the  factor  to  a 
solid  surface,  CONFAC  II  will  compute  the  factor  to  arbitrary 
solid  surfaces  or  regular  solids  such  as  parallepipeds,  cylinders, 
cones,  etc,,  with  the  restriction  that  all  of  the  surface  must 
appear  above  the  horizon  of  Surface  1,  CONFAC  I  cannot,  in  general, 
be  used  to  compute  factors  to  surfaces  which  are  occluded  or 
"shadowed"  in  a  varying  manner  by  intervening  surfaces;  on  the  other 
hand,  the  factor  in  such  instances  can  be  determined  by  CONFAC  II 
with  few  restrictions. 

6.  CONFAC  I  has  only  two  principal  classes  of  data  —  surface  data  and 
surface  transformation  data.  No  distinction  of  nsta  entry  is  made 
betvieen  plane  and  nonplane  surfaces.  Surface  data  is  distinguished 
from  transformation  data  by  the  position  of  the  data  name  on  tho 
data  name  card,  CONFAC  II,  however,  utilizes  nine  data  classifica¬ 
tions,  as  follows; 
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Class  1  -  Plane  polygon,  silhouette  developed  directly  from  data 

Class  2  -  Nonplane  polygon,  silhouette  developed  directly  from  data 

Ciass  3  -  Internally  generated  disk,  silhouette  developed  directly 
from  generated  data 

Class  4  -  Plane  polygon,  silhouette  internally  developed 

Class  5  -  Ncnpiane  surface  or  solid,  silhouette  ir*teisial]y  developed 

Class  6  -  Internally  generated  regular  disk  or  solid,  silhouette 
internally  developed 

Class  ‘7  -  Sphere 

Class  8  -  Multisurface,  silhouette  internally  developed  from  all 
surfaces  taken  together 

Class  9  -  Transformation  data 

7,  CON FAC  II  incorporates  a  silhouette  generator  subroutine  which  is 
utilized  when  the  factor  to  solids  or,  in  certain  cases,  to  non¬ 
solids  is  requested.  The  silhouette  generator  computes  the  per¬ 
spective  of  Surface(s)  2  from  preselected  positions  on  Surface  1 
from  which  configuration  factors  are  computed* 

8,  CONFAC  II  incorporates  an  internal  automatic  surface  generator 
which  computes  the  surface  boundary  coordinates  of  regular  plane 
and  solid  surfaces  from  input  data  specifications*  This  feature 
enables  the  analyst  to  create  surfaces  such  as  circular  or  ellipti¬ 
cal  disks,  parallelepipeds,  pyramids,  cones,  truncated  cones,  cylin¬ 
ders,  etc.  An  endless  variety  of  regular  surfaces  can  be  created 
by  CONFAC  II, 

9,  CONFAC  II  incorporates  extremely  fast  computation  of  factors  to  a 
sphere  which  is  arbitrarily  oriented  with  respect  to  Surface  1. 
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SECTION  II 


ANALYTICAL  PROCEDURES 


CONFIGURATION  AND  FORM  FACTOR 


The  general  eouation  that  must  be  solved  in  the  determination  cf  the 
radiant-interchange  fora  factor  is  (see  Figure  I) 


/ 
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cc  COS  01  COS  69  dA9dAi 


(1) 


The  following  part  of  the  integrand  is  the  factor  from  the  elemental  surface 
dA^  to  the  total  surface  A2,  referred  to  as  the  configuration  factor  or 
plane  point  factor,  e.p. 


c 
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cos  cos  02 


dA, 


(2) 


Therefore, 


(3) 


A  very  simple  geometric  interpretation  of  Equation  2  is  given  by  Nusselt 
The  principal  value  of  the  Nusselt  concept  is  that  the  computational  proce¬ 
dure  is  simplified  and  made  more  accurate  by  the  fact  that  no  mathematical 
or  numerical  integration  is  required  to  compute  the  configuration  factor.  How 
ever,  the  Nusselt  method  yields  only  the  configuration  factor  from  the  ele¬ 
mental  area  dA^ ;  one  must  still  integrate  all  such  factors  over  surface  Aj 
to  yield  the  form  factor  /^2  as  given  in  Equation  3, 

The  Nusselt  concept  utilizes  a  hemisphere  of  radius  R  constructed  over 
the  incremental  plane  area  dAj,  as  shown  in  Figure  1.  Every  point  defining 
the  boundary  of  surface  A?  is  projected  radiallv  to  the  hemisphere  surface 
and  then  vertically  downward  to  the  plane  of  dA^.  the  equatoiial  plane  of 
the  hemisphere.  The  locus  of  all  points  thus  projected  encloses  an  area, 

A‘*2>  CfTt  **Lc  hemisphere  haseA  This  area  A"2,  divided  by  the  area  of  the 
base,  is  the  configuration  factor  from  dAj  to  A2* 


The  validity  of  this  conclusion  can  be  demonstrated  as  follows.  Note 
that  the  elemental  area  dA2  is  described  in  surface  A-  by  the  elemental 
solid  angle  or 


dw 


1 


cos  02 

— 


(4) 


Similarly,  on  the  sphere  having  radius  R, 

dA2 

=  F* 


(5) 


Because  dA2  is  the  projection  of  dA2  on  the  hemisphere  base, 


,  dA2 

dAo  =  -  - 


2  ~  cos 

Inserting  Equation  6  in  Equation  5, 

tt 

d  Ao 


dw,  =  — - 

Rzcos  0 


(6) 


C7J 
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The  right  side  of  Equation. 4  appears  explicitly  in  Equation  1  and,  because 
Equation  7  is  identical  to  Equation  4,  Equation  2  becomes 


e  f  f  cos  0,  I  dA2  \ 

12  =  = 


For  a  sphere  of  unit  radius  (unit  sphere), 


JX2d#*  »; 


e  .JH. 

12  ir 


(8) 


which  completes  the  proof  of  Nusselt's  method.  By  inserting  Equation  8  in 
Equation  3„  t;»s  original  equation  becomes  greatly  simplified*,  only  one  area 
integration  is  now  required, 

>» 
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1  \  f  A2 

=  ^JJ  T-  dAl 


(9) 
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The  computer  program  described  herein  solves  Equation  9  numerically 
by  successive  algebraic  evaluation  of  A2  at  preselected  points  on  Surface 
Aj,  with  subsequent  numerical  integration  to  yield  /1?,  or 


f  1  TV  A? 

yi2  •  — 4^Al  T*  6A1  uo) 

•f 

It  should  be  emphasised  that  area  A2  is,  in  fact,  forced  by  the  doubly 
projected  silhouette  of  surface  A2  as  it  appears  from  dA^* 

The  element  dA^  is  assumed  to  be  oriented  in  the  XY  jjlane  and  at  the 
origin  of  the  coordinate  system  of  Surface  A2*  The  area  A?  can  be  found  „ 
from  the  line  integral  where  y^  =  F(x^)  is  the  locus  of  the  boundary  of  A2, 


A 


"  1 
2 *  ~r 


(Xi  dyj  -  yi  dxj) 


(11) 


Let  2  =  F(x,y)  be  the  locus  of  the  silhouette  of  A2>  and  S  the  distance  from 
dAj  to  the  point  (x,y,z)  on  the  silhouette  of  A2. 


S  =\l  x2  +  y2  +  z2 


Trcm  similar  triangles, 

X1  =  -7T'  <lxl  =  -5“  <1X  +  X  d(  '5— ) 

yl  =  ”1 T»  =  "T“  dy  +  y  d("T") 
Inserting  in  Equation  11 


*2  -4-/r 


xdy  -  ydx 
S2 


(12) 


Equation  12  can  be  transposed  to  finite  difference  form  by  replacing 
the  differentials  with  increments  for  numerical  evaluation.  Because  of  the 
problems  of  increment  size  control,  it  appeais  desirable  to  solve  Equation 
12  for  a  finite  line  segment  in  space  and  to  allow  the  analyse  to  control 
accuracy  of  configuration  factor  computation  by  suitable  selection  of  line 
segments  describing  Surface  2f  If  the  surface  is  actually  a  polygon  or  poly- 
hedra,  the  simulation  is  perfect;  if  the  surface  boundary  is  curved,  like  a 
disk,  for  example,  the  validity  of  the  result  is  a  function  of  the  number 
of  line  segments  used. 
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However,  a  much  simpler  and  more  easily  understood  geometric  d.  riva- 
tion,  using  the  unit  sphere,  yields  the  result  in  superior  computational 
form*  Referring  to  Figure  2,  note  that  the  radial  projection  of  line  seg¬ 
ment  AB  on  the  hemisphere  surface  forms  the  circular  arc  A:Bf.  Projection 
of  to  the  base  plane  produces  the  elliptical  arc  A,,B,,>  forming  the 
elliptical  section  A' '  0  B 1 1  with  the  origin* 

If  all  line  segments  describing  Surface  2  are  similarly  projected, 
the  area  ^  will  be  foimed  by  a  closed  series  of  elliptical  arcs.  Surface 
Az  does  not  have  to  he  a  plan**  Actually,  the  area  A2  results  from  the 
geometry  of  a  silhouette;  any  surface  or  object  projecting  an  identical 
silhouette  in  the  same  spatial  position  on  the  hemisphere  surface  will  pro¬ 
duce  the  same  area  A2  and  the  same  point  factor* 

11 

*nspection  reveals  that  the  magnitude  of  area  V>  can  ce  determined 
by  com),  ing  the  area  of  each  elliptical  sector,  properly  signed,  followed 
by  an  algv.  raic  summation. 

In  Figure  2,  the  area  of  elliptical  sector  Ae  is  the  projected  area 
of  circular  sector  Ag.  If  the  angle  between  the  plane  of  the  circular 
sector  A'OB*  and  the  XY  plane  is  y,  then 


The  area  Ag  is  computed  from  the  usual  polar  equation,  with  0  in  radians, 

1  2 
As  «  ~r  re 


For  the  unit  radius  sphere, 
As  ■  2 


(14) 


Substituting  Equation  14  in  Equation  13,  and  solving  for  Ae, 


Ac  =  —5—  cos  y 


(15) 


For  a  polygoT.  uf  N  sides,  the  net  area  A2  is  found  by  algebraic  summation 
of  all  computed  A£. 
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A2  =  -y- 


CCS  Y 

n  ’n 


n=l 


(16) 


Substituting  in  Equation  8,  we  have 

i  N 


*12  “  ~k  j  2  6n  C0S  Y"  | 

I  n=l  1 


(17) 


A  general  analytical  derivation  of  this  equation  is  given  in  Reference 
3,  and  is  reported  to  be  originally  developed  by  Omoto  in  1924. 

The  absolute  value  notation  will  be  explained  later.  The  use  of  vector 
algebra  greatly  facilitates  the  consultation^  of  0  and  cos  y*  Taking,  for 
example,  directed  line  segments  of  OA  and  OB,  the  vector  dot  product  is 


OA  •  OB  =  xAxB  +  yAyB  +  zAzB 


(18) 


The  cross  product  OA  x  OB  in  determinant  form  is 

i  j  k 


OA  x  OB 


*A  '/A  ZA 
xb  Yn  ZB 


which,  upon  expansion,  becomes  the  normal  vector  V^., 

VN  =  OA  x  OB  =  (yAzB  -  zAyB)i  +  (xBzA  -  zBxA)j  +  (xAyB  -  xByA)k  (19) 


where  i,  j,  and  k  are  mutually  orthogonal  unit  base  vectors  directed  along 
the  principal  axes. 

is  equal  in  magnitude  to  twice  the  area  of  the  triangle  AOB  and 
is  oriented  normal  to  the  plane  of  AOB  so  that  the  three  vectr^s  form  a 
right-handed  system.  The  magnitude  is  computed  by  the  Pythagorean  theorem, 
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\%\  =  V 


(Vb  -  zAyB>  +  <x 


BZA 


Vb1  ' 


'Vb  -  xByA^  ‘ 


(20) 


The  angle  8  may  be  evaluated  from  either  the  dot  or  the  cross  product  by 
use  of  inverse  functions,  specifically 


0 


OA  *  OP 

to  a' I  fi® 


.  -1 
or  sin 


ToaJ  !ob| 


1 

! 


‘however,  an  overall  economy  of  computation  results  from  the  use  of  the 
arctan  function. 


0  =  tan 


-1 


VN1 


OA  •  OB 


(21) 


As  noted  earlier,  the  angle  y  is  defined  as  the  angle  between  the 
plane  of  AOB  and  the  XY  plane.  It  is  also  the  angle  between^tne  vector 
and  the  2  axis;  cos  y  is  therefore  the  direction  cosine  of  with  respect 
to  the  2  axis.  Using  the  Z  component  in  Equation  19. 


cos  y 


xAyB  *  xByA 


(22) 


If  the  numerator  and  denominator  are  both  divided  by  2, 


cos  Y  = 


WB  ^  XBYA 

m 


This  shows  that  cos  y  is  also  equal  to  the  ratio  of  the  signed  projected  area 
of  triangle  AOB  on  the  XY  plane  and  the  plane  area  of  triangle  AOB. 


In  the  right-handed  system  shown,  cos  Y  is  positive  when  the  order  of 
computation  of  the  vectors  in  the  cross  product  causes  the  normal  vector 
to  point  in  the  direction  of  the  +Z  axis  (O<y<90),  The  order  in  which  one 
proceeds  from  point  to  point  on  the  boundary  of  Surface  2  wil t  sign  each 
elliptical  sector  accordingly;  however,  because  the  sector^  are  summed 
algebraically,  the  same  absolute  magnitude  will  result  regardless  of  order* 
Because  the  point  factor  is  always  a  positive  number,  the  order  is  computa¬ 
tionally  unimportant.  Nevertheless,  the  program  requires  that  data  be 
entered  in  counterclockwise  order  for  other  reasons.  This  will  be  dis¬ 
cussed  in  more  detail  later. 
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The  relative  ease  with  which  the  point  factor  can  be  computed  is  best 
illustrated  by  an  example.  Using  the  triangle  shown  in  figure  2,  and  start¬ 
ing  with  line  segment  AB,  from  Equation  18 

3A#0B=l+3+9=13 


from  Equation  20 


I  =  I«A  x  OB  j  =^-6)2  +  0  +  (2) 2  =i/40“ 


From  Equation  21 


From  Equation  22 


8AB  =  tan'1  [^]  ?  0.453 


'  cos  yab  %75r  =  °*316 


Moving  to  BC, 

OB  •  OC  =  3  +  3  +  9  =  15 


cos  ynr  =  =  -0,686 

L  VI36 


13 


Finally,  line  segment  CA 


0C*0A  =  3+1+9  =  13 
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CA 


0*453 


cos 


0.316 


The  configuration  factor  is,  therefore,  from  Equation  17, 


C12  =  “ZiT  j  2(0.453)  (0.316)  +  (0.661)  (-0.686) 

=  T7  |  -°-167| 

C12  -  0.0266 


Note  the  repetitive  nature  of  the  computation.  Thus,  all  surfaces 
represented  by  straight  line  segments  in  space  can  be  analyzed  in  the  simple, 
direct  manner  shown. 


COORDINATE  TRANSFORMATION 

The  task  of  computing  factors,  even  when  simple  "closed- form"  solutions 
are  available,  is  oftentimes  laborious  because  the  surfaces  under  considera¬ 
tion  appear  in  difficult,  skewed  relative  positions.  A  significant  nart  of 
this  effort  has  been  eliminated  by  the  program  through  the  canabilitv  of  gen¬ 
eral  coordinate  transformation  (translation  and/or  rotation)..  Surface  data 
may  he  entered  for  each  surface  using  an  individually  convenient  local  origin. 
The  surfaces  may  then  be  linked  together  by  transforming  one  or  both  surfaces 
to  a  convenient  third  origin  which  is  common  to  both  surface*.  The  fact  that 
internally  generated  surfaces  may  edso  be  transformed  (excluding  multisurfaces; 
makes  this  feature  a  very  powerful  tool. 
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Actually,  two  different  types  o f  coordinate  trsnsforinritic 
the  program.  The  transfonnation  discussed  in  the  prior  paragraph  is  termed 
a  "primary11  transfonnation,  and  is  under  control  of  the  user  through  trans¬ 
formation  data  entry.  The  second  type  of  transformation  is  termed  ar.  "auxil¬ 
iary"  transformation, and  is  unoer  internal  program  control  only.  An  auxiliary 
transformation  transforms  the  surface  coordinates  of  both  surfaces  into  a  new 
coordinate  system  formed  so  that  the  XY  plane  of  the  coordinate  system  lies  in 
the  reference  plane  of  one  of  the  surfaces.  The  reference  plane  of  a  surface 
is  the  plane  formed  by  the  first,  second  and  last  point  describing  that  sur- 
fa  e.  The  origin  of  an  auxiliary  coordinate  system  is  located  at  point  i  in 
the  particular  surface  controlling  the  transformation.  The  X-axis  is  direct¬ 
ed  along  the  line  segment  formed  by  points  1  and  2,  The  surface  unit  orienta¬ 
tion  vector  becomes  the  Z  axis;  the  Y  axis  is  computed  orthogonal  to  the  X  and 
Z  axes,  thus  locating  the  XY  plane  in  the  control  surface  reference  plane. 


The  auxiliary  transformation  actually  serves  two  purposes.  First,  it  is 
utilized  by  Subroutine  DOICU  to  facilitate  reconstruction  of  the  "seen"  part 
of  surfaces  which  are  not  entirely  seen  by  the  other  surface.  Secondly,  the 
program  requires  that  prior  to  computation  of  the  configuration  factors,  Sur¬ 
face  1  must  appear  in  the  XY  plane  of  the  final  coordinate  system  alone  with 
Surface  2  in  its  proper  relative  position.  This  is  necessary  to  enable  Sub¬ 
routine  MAP  to  select  points  on  Surface  1  from  which  factors  to  Surface  2  may 
be  directly  computed,  or  from  which  silhouettes  of  Surface  2  may  be  generated 
and  factors  computed. 


For  example,  suppose  Figure  3  represents  the  surfaces  of  various  items  of 
equipment  appearing  in  a  compartment.  The  unprimed  coordinate  system  shown  may 
be  conveniently  chosen  at  a  comer  or  axis  of  symmetry,  perhaps  as  shown  on  a 
mechanical  drawing,  This  system  may  not  be  convenient  for  data  entry  of  the 
disk,  however.  The  primed  coordinate  system  with  the  origin  at  the  center  of 
the  disk  is  the  more  logical  choice  in  this  case.  The  previously  described 
surface  generator  will  generate  the  disk  about  this  origin.  The  disk  data  can 
then  be  transformed  from  the  primed  to  the  unprimed  system  by  a  primary  trans¬ 
formation,  The  choice  of  generating  the  cube  and  transforming,  or  directly 
entering  data  from  the  unprimed  system,  is  left  to  the  user  as  it  requires 
about  equal  effort  both  ways.  The  plate  coordinates  can  be  easily  entered 
from  the  unprimed  system.  Now,  suppose  we  desire  the  form  factor  from  the 
disk  to  the  plate.  If  the  data  are  entered  as  discussed  above  (including  the 
transformation  data),  the  program  will  generate  the  disk  and  then  primary 
transform  disk  coordinates  to  the  unprimed  system.  Since  the  disk  is  bisected 
by  the  plate,  an  auxiliary  transformation  of  all  coordinates,  both  disk  and 
plate,  will  be  made  from  the  unprimed  to  the  quad-primed  system.  Now,  that 
portion  of  the  disk  appearing  above  the  active  side  of  the  plate  wilt  be  de¬ 
termined,  and  an  auxiliary  transformation  of  the  plate  and  the  truncated  disk 
will  bo  made  to  the  double  primed  coordinate  system,  i.e,,  the  reference  plane 
of  the  disk.  The  disk  is  now  in  a  position  for  mapping,  and  the  plate  coor.« 
dinates  are  proper  for  obtaining  the  configuration  factors,  A  similar  manip¬ 
ulation  of  surface  data  would  be  made  to  obtain  the  form  factor  to  the  cube 
with  one  exception  -no  truncation  of  the  disk  would  occur  and  the  auxiliary 
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transformation  lo  the  double -primed  system  would  occur  iuITiCdlutC  ly  after  uri¬ 
nary  transformation* 

The  transformation  technique  utilized  for  a  primary  transformation  dif¬ 
fers  from  the  customary  method  whereby  ”old"  coordinates  plus  translation 
data  and  direction  cosines  or  Euler  angles  are  supplied,  from  which  a  "new1* 
set  of  coordinates  are  derived*  The  program  requires  the  coordinates  of  any 
three  points  (not  in  a  line)  measured  from  the  new  origin.  These  data  arc  then 
used  to  derive  direction  cosines  and  translation  terms,  hy  which  the  old  coor- 
dir.aces  are  then  transformed  to  the  new  origin. 

The  reader  may  find  it  easier  to  visualize  transformation  in  terms  of  the 
movement  of  the  surfaces  instead  of  the  origins.  In  the  case  of  the  disk, 
again  referring  to  Figure  3,  we  may  say  we  generated  the  disk  with  its  center 
at  the  origin  of  the  unprimed  system  and  in  its  XY  plane,  .and  then  moved  the 
surface  to  the  position  indicated  by  the  primed  system.  This  viewpoint  appears 
more  realistic  when  motion  is  simulated  by  transforming  a  surface  along  a  par¬ 
ticular  path. 

The  mathematical  treatment  of  primary  and  auxiliary  transformation  is 
presented  in  Appendix  C. 


SILHOUETTE  GENERATOR 

As  noted  in  the  introduction,  CONFAC  I  cannot,  in  general,  be  used  to 
compute  'the  form  factor  to  solid  surfaces.  Subroutine  FACTOR  requires  a  sin¬ 
gle  array  be  made  available  containing  the  surface  boundary  points,  and  only 
those  points,  which,  when  taken  in  numerical  sequence,  form  a  valid  silhouette 
of  Surface  2  from  a  particular  point  in  Surface  1,  It  is  impossible  to  gen¬ 
erally  satisfy  this  requirement  with  a  single  input  array  if  Surface  1  is  finite 
and  Surface  2  is  arbitrarily  nonplanar  or  solid.  It  is  the  function  of  the 
silhouette  generator  to  determine  which  points  in  a  given  set  cf  Surface  2  data 
form  the  silhouette  fiom  preselected  viewpoints  on  Surface  1. 

The  silhouette  generator  computes  the  silhouette  from  the  perspective  of 
Surface  2  developed  on  the  Z-unity  (Z  »  +1)  plane.  The  perspective  on  the  Z- 
unity  plane  is  the  locus  formed  in  the  plane  by  the  boundary  of  the  solid  angle 
subtended  by  Surface  2, 

For  example,  the  view  of  a  cube  from  two  positions  on  tho  XY  plane  is 
shown  in  Figure  4,  The  coordinates  of  each  point  in  the  Z-unity  plane  are 
derived  in  the  following  manner  from  the  coordinates  of  its  corresponding 
point  on  the  cube.  Note  the  triangle  formed  by  the  origin,  point  2  in  the  cabe, 
and  point  Q,  the  vertical  projection  of  point  2  on  the  XY  plane,  A  similar 
triangle  is  constructed  from  point  2*  to  point  N,  From  similar  triangles, 

ON  Z2 

m  *  ?2  (22) 
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TOP  VIEW 


FIGURE  4.  SIMPLE  SILHOUETTE  GEOMETRY 


In  like  manner,  using  triangles  RON  and  SOQ  in  the  ton  view, 


but  OR  =  X£,  OS 
Therefore, 


ON  OR  _  NR 

00  =  0?  =  Q? 

X,,  NR  =  Y2<  QS  =  Y2  and  z\  =  1. 
♦  ♦ 

x2  Z2  _  1 

x2  =  27 


h 

“2 


(23) 


and  similarly, 


(24) 


This  reduction  to  two  dimensions  results  in  considerable  simplification. 
Given  the  coordinates  (X,  Y)  and  point  connections  data,  it  is  possible  to 
determine  the  line  segments  forming  the  silhouette  by  application  of  a  simple 
criterion,  At  each  point  on  the  silhouette,  those  line  segments  forming  the 
largest  included  angle  define  the  silhouette.  For  example,  at  point  2f  in 
the  lower  silhouette  m  Figure  4,  vctLuxs  4.^  -  o  ,  4.  -1  <uiu  c  -  4  emerge 

from  the  point.  Vectors  2  -  8f  and  2*  -4*  obviously  form  the  silhouette, 

and  can  be  numerically  selected  by  applying  the  criterion. 


Figure  5  shows  the  development  of  the  Z-unity  plane  silhouette  of  a 
multisurface  .  In  contrast  to  point  D,  Surfaces  S2  and  S3  apnear  separated 
in  the  silhouette  when  viewed  from  C. 


Note  the  line  connecting  4  to  7  ,  This  artifice  -a  "bridge*1  line  -  is 
utilized  to  cause  the  silhouette  generator  to  include  both  surfaces  in  the 
silhouette,  otherwise  surface  S3  would  he  ignored.  Because  the  line  has  no 
width,  it  has  no  effect  on  the  factor  computation,  but  the  silhouette  gener¬ 
ator  follows  the  line  as  if  it  were  a  boundary  of  the  multisurface  S2  plus  S3. 

The  distinguishing  difference  between  the  silhouettes  shown  m  Figure  4 
and  Figure  5  is  the  fact  that  "crossover"  occurs  in  Figure  5.  The  silhouette 
at  a  crossover  is  formed  by  intersecting  line  segments  at  a  point  between  line 
segment  extremities.  The  detection  of  such  intersections,  and  the  computation 
of  the  coordinates  of  the  intersection,  requires  considerable  analysis  with 
resultant  increased  computer  time,  Rccauso  of  this,  siJhouetVv*  analysis  is 
termed  "simple"  if  no  investigation  is  made  by  the  silhouette  generator  to 
detect  crossovers,  and  "complex"  when  such  is  made.  Only  multisurface  data 
(class  8)  are  run  ill  the  complex  mode.  All  other  surface  data  requiring  the 
silhouette  generator  (classes  4,  5  and  6)  are  run  in  the  simple  mode. 


SURFACE  GENERATOR 


The  urogram  cannot directly  compute  factors  to  curved  surfaces  or  bound¬ 
aries  such  as  disks,  cylinders,  etc.  A  series  of  line  segments  must  be  sub¬ 
stituted  for  a  curved  line*  In  general,  the  more  line  segments  used,  the  more 
accurate  the  simulation*  Because  every  surface  ooint  reouires  3  coordinates 
(and  connecting  joints  data,  when  the  silhouette  must  be  computed),  prepara¬ 
tion  and  entry  of  data  for  even  a  modest  simulation  of  a  cylinder  can  involve 
a  considerable  amount  of  effort*  The  internal  surface  generator  eliminates 
practically  all  of  this  effort. 

The  surface  generator  is  used  to  create  surfaces  entered  under  data  clas¬ 
ses  3  or  6*  Regular  plane  polygons  are  created  under  Class  3,  but  no  connec¬ 
tions  data  are  generated.  A  regular  plane  polygon  or  solid  surface,  including 
connections  data,  is  created  under  a  class  6  entry. 

The  surface  generator  "creates"  a  surface  in  accordance  with  cross  section 
specifications.  The  following  information  is  required  to  create  a  class  3 
surface: 


1.  Number  of  cross  section  division  (sides)  >  3 

2.  Coordinates  (X,  Y,  Z)  of  center  of  polygon 

3.  X-axis  Radius 

4*  Y-axis  Radius 

Because  a  class  6  surface  may  have  one  or  more  cross  sections,  the  following 
data  are  required: 

1,  Number  of  cross  section  division  (sides)  >  3 

2,  Number  of  cross  sections 

3,  Coordinates  (X,  Y,  Z)  of  first  cross  section 

4,  X-axis  Radius  of  first  cross  section 

5,  Y-axis  Radius  of  first  cross  section 

If  more  than  one  cross  section  is  specified,  the  following  data  are  renuired 
for  each  additional  cross  section:  X-axis  radius,  Y-axis  radius  and  Z-cocr- 
dinate.  All  cross  sections  are  created  parallel  to  the  XY  pla^o  of  the  gen¬ 
erator  coordinate  system,  and  must  be  specified  above  the  XY  plane*  Note  that 
X,  Y  coordinates  are  required  to  locate  "he  first  cross  section  only.  If 
more  than  one  is  specified,  all  are  oriented  along  the  same  vertical  center- 
line  to  the  position  specified  by  the  respective  Z  coordinate. 


The  basic  generating  element  is  the  ellipse.  Because  only  complete 
polygons  are  generated,  the  total  angle  of  2tt  radians  about  the  vertical  cen¬ 
terline  is  divided  by  the  number  of  sides  snecified  to  yield  the  unit  parametric 
angle  $  in  the  equations  of  the  ellinse: 


N 

X  =*  (XR)  cos  $ 
Y  =  (YR)  sit.  <J> 


Instead  of  the  conventional  semi-major  and  semi-minor  expressions,  the 
terms  "X-Radius*’  and  "Y-Radius"  are  utilized  -  the  larger  of  the  two  becomes 
the  semi-major  axis  as  shown  in  Figures  6  (a),  (b) ,  and  (c) ,  When  XR  -  YP, 
the  generating  figure  is  circular  and  a  regular  polygon  of  N  sides  results. 
Notice  that  the  generating  figure  always  circumscribes  the  generated  polygon. 
The  radius  vector  always  starts  in  the  same  relative  position  parallel  to  X- 
axis  and  moves  counterclockwise  about  the  vertical  centerline  of  the  generated 
cross  section.  Considerable  computing  time  is  saved  by  using  sin  (4>  +  S),  cos 
(4>  +  8)  trigonometric  formulae  for  computation  of  X,  Y  after  unit  values  are 
obtained  by  use  of  computer  library  functions. 

Figure  6  (d)  shows  an  eight-sided  polygon  elevated  above  the  XY  plane. 
Figure  6  (e)  indicates  the  order  in  which  point  numbers  are  assigned  to  a 
solid  surface.  The  first  point  is  always  assigned  to  the  first  coordinates 
in  the  first  cross  section.  Numbers  are  assigned  in  numerical  sequence  ver¬ 
tically  until  the  last  cross  section  is  numbered,  for  a  particular  value  of 
the  sequence  is  continued  in  similar  manner  with  the  first  cross  section  and 
the  next  value  of  until  all  points  are  defined. 

Point  connections  data  are  also  computed  for  each  point  for  Class  6  sur¬ 
faces,  For  example,  in  Figure  6  (c) ,  points  2,  3  and  15  are  computed  for 
point  1;  points  1,  4  and  16  are  computed  for  point  2,  etc,  This  information 
is  used  by  the  silhouette  generator. 

An  example  of  the  variety  of  objects  which  may  be  created  by  a  few  cards 
of  specifications  are  shown  in  Figure  6  (f) ,  The  cone  vertex  is  generated 
merely  by  specifying  zero  X-radius  and  Y-radius, 

The  internal  surface  generator  also  commutes  the  surface  area  of  the  gen¬ 
erated  scl-5 d 5  if  the  cross  sections  are  similar.  Because  thr  silhouette  gen¬ 
erator  analyzes  the  solid  figure,  the  total  surface  are:,  is  cemented.  For 
instance,  the  area  olHtKe  prismatic  cylinder  shown  in  Figure  6  (cj  would 
include  the  top  and  bottom  polygons.  The  surfaco  area  computation  analytical 
development  is  given  in  Appendix  D. 
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SECTION  III 


COMPUTER  PROGRAM  CONTENTS 


PROGRAM  DESCRIPTION 

The  program  is  written  in  iBM  7090  FORTRAN  II  source  language.  The  source 
deck  consists  of  the  Main  Program  and  Subroutines  UNIVEC,  TXFRM,  COICU,  MAP, 
SILFAC  and  FACTOR,  An  input-output  tape  compatibility  Subroutine  written  in 
IBM  7090  (FAP)  machine  language  is  included.  Algebraic  routines  required  from 
lit ‘-arv  tape  are  SORT  (Square  Root),  ARC7AN  (Inverse  tangent  trigonometric 
function),  COS  (Cosine  trigonometric  function),  and  SIN  (Sine  trigonometric 
function).  The  source  programs  are  presently  dimensioned  so  that  a  32  K  core 
size  is  required,  NAA  Library  Subroutines  COUNTV  and  T1MEV  are  also  used  by 
Subroutine  SILFAC  when  operating  in  the  NAA  7094  system. 

Main  Program 

The  functions  performed  by  the  Main  Program  are  as  follows: 

1,  Reads  in  surface,  transformation  and  run  data, 

2,  Processes  innut  surface  data  and  prints  immediately  unon  completion. 
Run  instruction  data  are  read  in  and  processed  one  card  at  a  time 
and  processed  at  once.  No  printout  of  the  complete  run  instructions 
is  ‘given,  as  was  in  COiTFAC  I. 

3,  Selects  the  proper  data  for  processing  according  to  the  run  ir>$truc 

t i one  ^ 

4,  Examines  each  run  instruction  and  calls  in  proper  subrout ines  for 
processing, 

5,  Prints  diagnostic  error  indications  when  possible, 

6,  Prints  standard  or  detailed  output  as  indicated  by  run  instructions. 
Subroutine  UNIVEC 


This  subroutine  computes  the  components  of  a  unit  orientation  vector 
normal  to  the  reference  plane  formed  by  the  first,  second  and  last  point  in 
surface  data  classes  1,  3,  4  and  plane  6,  The  cross  product  of  vectors  1-2 
and  1-last  is  computed  and  normalized.  The  vector  is  formed  normal  to  point 
1,  and  is  located  on  the  active  side  of  the  surface,  thus  orienting  the  sur¬ 
face. 

It  also  computes  a  new  fourth  point  normal  to  the  new  three  points  sub¬ 
mitted  in  transformation  data  and  an  old  fourth  point  normal  to  the  old 
three  points  in  the  surface  data  to  be  transformed. 
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Subroutine  TXFRM 


The  first  section  perforins  the  auxiliary  transformation*  This  trans¬ 
formation  is  used  to  reconstruct  a  surface  which  is  bisected  by  the  second 
surface*  It  also  tests  Surface  i  to  determine  if  the  reference  plane  is  sub¬ 
stantially  in  the  XV  plane  of  its  coordinate  system.  If  it  is  not,  an  aux** 
iliarv'  transformation  is  effected  to  move  the  surfaces  to  fulfill  this  re¬ 
quirement  nrior  to  computation  of  silhouettes  or  factors* 

This  subroutine  also  nerfoms  a  primary  transformation  as  indicated  by 
?n.  instructions  and  t ran s format ion  data*  This  transformation,  if  indicated 
for  a  surface,  is  accomplished  nrior  to  entry  to  subroutine  DOICU  so  that  tests 
of  the  surface  "view"  of  each  other  occur  in  their  transformed  position (s). 

Subroutine  DOICU 


The  function  of  this  subroutine  is  conveyed  literally  by  its  name  DO-I- 
C-U*  (liven  surfaces  A1  and  A2  with  the  "active"  side  of  each  surface  identi¬ 
fied  by  the  surface  orientation  unit  vector j  the  question  is  asked;  Is  all, 
part,  or  none  cf  surface  A1  "seen"  by  A2?  Conversely,  does  A2  see  all,  none, 
or  part  of  Al?  This  is  accomplished  by  computing  the  vector  dot  product  form¬ 
ed  by  the  unit  vector  in  one  surface  with  the  vector  formed  by  point  1  in  the 
first  surface  and  each  point  in  the  other  surface  (See  Figure  7).  The  sign  of 
the  dot  product  indicates  whether  the  angle  between  the  vectors  is  less  than 
or  greater  than  90° ,  which  reveals  the  position  of  the  point  relative  to  the 
of  the  viewing  surface.  In  Figure  7  (a)  the  dot  products  from  surface 
Al  to  A2  are  all  positive,  and  conversely,  all  from  A2  to  Al  arc  likewise  pos¬ 
itive;  Al  sees  all  of  A2;  A2  sees  all  of  Al,  However,  in  Figure  7  (b)  all 
dot  products  from  A2  to  Al  are  positive,  but  all  from  Al  to  A2  are  negative. 
Hence,  in  general,  if  all  dot  products  from  one  surface  to  another  are  all 
negative,  then  the  surfaces  do  not  see  each  other,  even  though  the  converse 
products  may  be  positive.  There  is  also  the  trivial  case  where  all  products 
are  zero,  in  which  case  the  surfaces  are  in  the  same  plane,  and  obviously 
again  cannot  see  each  other. 

Figure  7  (c)  shows  a  surface  A2  bisecting  surface  Al,  In  this  case, 
some  of  the  dot  products  from  A2  and  Al  are  positive  and  some  negative.  In 
Figure  7  (d)  both  Al  and  A2  are  bisected,  Nonnlanar  surface  A3  was  added  to 
show  how  it  would  be  bisected  by  Al*  Surface  A3  has  no  orientation  vector  and 
thus  no  test  is  made  of  the  view  from  this  surface.  The  vertical  dashed  line 
in  A2  represents  how  the  plane  1-2-5-6  in  A3  might  bisect  A2,  DOICU  will  not 
detect  this  condition.  If  the  configuration  factor,  C ^  were  required,  DOICU 
would  properly  bisect  A3,  However,  if  the  factor  to  the  concave  side  oplv  is 
desired,  an  error  would  result  because  part  of  A2  sees  the  convex  side  of  A3. 
This  represents  one  of  the  limitations  of  CONFAC  I  which  canied  over  to 
CON FAC  ir. 

If  a  surface  is  bisected,  DOICU  reconstructs  the  surface  data  to  exclude 
the  area  not  scon  by  the  other  surface.  If  point  1  in  the  original  surface  is 
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removed  as  a  result,  a  new  orientation  vector  is  created  ovc^  the  new  point  1 
as  shown  in  Figure  7  (c).  Notice  that  in  reconstructing  A3  [Figure  7  (d)], 
DOICU  created  the  new  array  1,2,  3,  4,5,  This  ‘'surface'*  is  identical  to 
the  actual  surface  seen  by  A1  insofar  as  factor  computation  from  At  is  con¬ 
cerned. 

The  bisection  of  a  surface  is  done  in  a  simple  manner,  with  the  aid  of 
the  auxiliary  transformation  capability.  For  example,  in  Figure  7  (c),  the 
coordinates  of  both  surfaces  are  transformed  so  that  A2  lies  in  the  XY  plane 
of  the  auxiliary  (primed)  coordinate  system.  Each  point  in  A!  is  tested r  in 
numerical  order,  until  a  change  in  the  sign  of  the  Z-coouIIiiate  occurs.  The 
coordinates  of f the  new  points  where  the  transition  line  segment  crosses  the 
X  Y  plane  (Z  =  0)  are  obtained  by  computing  X  and  Y  intercepts  of  traces 
projected  on  X  7,\  Y1  Z1  principal  planes. 

Subroutine  MAP 


The  double  integral  in  Equation  (9)  and  its  numerical  counterpart  in 
Equation  (10)  mathematically  represent  the  volume  under  a  surface  defined  by 
the  configuration  factor  2  =  f(X,Y).  Subroutine  MAP  decides  the  location 
(X,Y)  from  which  each  factor  to  Surface  2  will  be  computed. 

It  is  assumed  that  Surface  1,  being  classed  as  a  plane,  is  a  plane  sur- 
face  throughout.  The  program  insures  only  that  the  reference  plane  of  Surface 
1  is  in  the  W  plane  of  the  final  coordinate  system.  MAP  will  use  the  X,Y 
coordinates  of  all  points,  and  assumes  a  value  of  0  for  all  Z  cooiuinates. 

This  procedure  cannot  properly  map  a  nonnlancr  surface. 

Subroutine  MAP  determines  the  maximum  Y  coordinate  and  the  minimum  Y  co¬ 
ordinate  from  among  the  points  defining  Surface  1  (Figure  8) .  The  total  ver¬ 
tical  distance  between  Y^ax  and  Y^^  is  divided  into  equal  vertical  increments 
as  specified  by  the  run  instructions.  Then,  horizontal  lines  are  scribed 
across  (parallel  to  X-axis)  the  surface  at  each  vertical  increment  position, 
including  Ymx  and  Yn^.  The  point  at  \t)  ich  horizontal  line  intersects  the 
left  (toward  the  negative  X  direction)  boundary  of  Surface  1  is  termed  "X- 
lcft"  and  the  intersection  on  the  right,  "X-right".  Each  horizontal  line  seg¬ 
ment  thus  created  is  termed  a  "mapping  line".  Each  mapping  line  segment  is 
also  divided  into  an  equal  number  of  increments  as  specified  by  the  run  instru 
tions.  All  mapping  lines  are  divided  into  the  same  number  of  increments,  not 
necessarily  the  same  size  of  increment.  Obviously,  if  Surface  1  converges  to 
n  point  instead  of  a  llmTat  Y  x  or  Y ,  the  horizontal  increment  is  0,  A 
configuration  factor  is  computed  at  each  increment  point  along  a  mapping  line, 
including  X-lcft  and  X-right,  which  means  the  number  of  factors  ner  line  is 
one  greater  than  the  number  of  increments. 

The  number  of  increments  is  automatically  set  to  24  horizontal  :nd  24 
vertical,  but  can  be  separately  specif:  by  input  data  to  3,  6,  18,  24,  30, 
36,  42,  48,  54,  or  60.  The  details  are  discussed  in  Section  IV, 
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A  typical  example  of  Surface  1  manning  using  a  standard  (24  x  24)  incre¬ 
ment  is  shown  m  Figure  8  (c),  The  manning  area  is  also  confuted  bv  Subroutine 
MAP,*  it  is  the  snm  of  the  rectangular  areas  formed  bv  each  manning  line,  A 
measure  of  form  factor  accuracy  is  the  degree  with  which  the  manning  area 
approximates  the  actual  surface  area. 

Figure  8  (a)  illustrates  a  Surface  1  orientation  which  cannot  be  satis¬ 
factorily  manned  because  the  crosshatched  area  is  ignored.  The  program  does 
not  detect  more  than  one  left  and  one  right  intersection  between  a  manping 
■  *  e  and  the  su'*iace  *-~>r'd3i  ,  therefore,  n"i«t  3  is  ignored.  "Hie  same  surface 
rotated  sufficiently  may  be  accentable,  however,  clearing  this  restriction,  as 
shown  in  Figure  8  (b)., 

S jb routine  SILF AC 


This  subroutine  computes  the  silhouette  of  Surface  2  which  annears  from 
the  points  selected  on  Surface  1  by  Subroutine  MAP,  and  then  computes  the 
configuration  factor  from  this  silhouette.  After  all  configuration  factors 
have  been  determined,  the  form  factor  is  comnuted  bv  numerical  integration. 

Surface  data  entered  as  Class  4,  5,  6,  7  or  8  is  processed  by  SILFAC. 

Class  7  data,  a  sphere,  is  processed  in  this  routine,  but  the  silhouette 
generator  is  not  utilized;  a  closed  form  solution  is  usea  instead  (see 
Appendix  E) . 

Classes  4,  S  and  6  are  processed  bv  the  silhouette  generator  in  the 
simple1  mode;*  only  those  points  given  in  connections  data  are  analyzed  to 
select  the  next  point  on  the  silhouette.  Surfaces  such  a*  planes,  cylinders, 
paralleluepids,  etc,  may  be  processed  in  the  simple  mode, 

Class  8,  a  multisurface,  is  the  only  class  processed  in  the  complex 
mode.  One  or  more  (limitations  on  data  are  given  at  the  end  of  this  section 
and  in  Section  IV)  surfaces  mav  be  processed  as  a  multisurface.  Processing 
m  the  complex  mode  is  complicated  because  the  computer  must  test  all  line 
segments  m  all  surfaces  (including  the  surface  in  which  the  segment  appears) 
in  order  to  select  the  next  point  forming  the  silhouette,  or  to  compute  the 
next  point  on  the  silhouette.  This  analysis  is  further  complicated  by  ambi¬ 
guities  insulting  from  normal  imprecision  in  input  data  and  internal  arithmetic 
roundoff,  necessitating  the  use  of  numerous  time  consuming  tolerance  tests. 
Consequently,  a  factor  computation  m  the  complex  mode  takes  considerably 
longer  than  the  simnlc  mode 

If  a  surface  is  processed  as  simple  when  it  should  be  comnle*  ,*  A 
wrong  silhouette  will  be  computed  whenever  a  crossover  (an  intersection  of 
two  lire  seg<icats)  occurs.  The  configuration  factor  confuted  from  that  sil- 
houc' to  will  be  wrong. 

It  is  possible  to  detect  certain  hinds  of  trouble  in  the  silhouette  gen¬ 
erator  from  the  detailed  silhouette  output  which  lists  the  points  forming  the 
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silhouette.  Normally,  the  silloiette  will  start  at  the  lowest  leftmost  roint 
in  the  Z-umtv  nlane  nersnectiv'*  devcloned  for  tne  noted  nomt  in  Surface  1, 
and  move  progressively  from  nomt  to  noint  m  a  counterclockwise  manner,  keen¬ 
ing  the  perspective  area  to  the  left.  When  a  crossover  occurs,  the  nomt 
is  computed  and  assigned  the  ne  t  Highest  number  in  the  silhouette  arrav.  For 
cxa*rplc,  from  point  "n"  m  Pipure  5,.  the  silhouette  derived  m  the  cnn-lcx 
node  would  annear  in  the  detailed  silhouette  outnut  as  "Line  No.,  Point  No,, 

1»  2',‘  O',*  6*,'  S’,’  S"  10','  47  1."  However,  if  this  nrohlen  were  run  as  simnle 
(both  surfaces  entered  together  as  Class  4)  instead.  Surface  S3  would  he 
i/no/cd  because  f*  crn«sso'"‘t*  at  noint  0  would  not  be  commuted.  The  sil¬ 
houette  would  annear  normal,  hut  actual lv  he  wrong,  as  follows 

Line  No.,  Point  No.,  1’,’  27  37  4','  lV 

A  bad  silhouette  can  sometmes  be  detected  bv  the  nresence  of  internal  "loop¬ 
ing".  Normally,  a  silhouette  is  completed  bv  a  return  to  the  starting  noint. 
^ut,  if,  for  some  reason,  a  wrong  nath  is  chosen,  it  nav  loon  a  nolcgon  within 
the  nersnectivc.  Loonint*  is  characterized  bv  the  rencated  appearance  the 
same  sequence  of  numbers.  No  internal  pattern  recognition  is  attempted,  the 
onlv  detection  i»  visual  examination  of  the  detailed  outnut. 

The  coordinates  of  the  silhouette  on  the  Z  =  l  nlane  are  used  directlv 
for  factor  commutation  instead  of  the  actual  noints  on  the  surface  m  space,. 
Because  the  Z  -  coordinate  of  each  noint  is  1^  the  configuration  factor  eoua- 
tions  for  this  snecial  case  can  be  simplified,  and  commuting  time  reduced, 

SI LT AC,  tlieiefuie,  contains  its  own  enualiuns  for  configuration  factor  commuta¬ 
tion  and  numerical  integration  across  Surface  1,  The  integration  nrocess  is 
similar  to  th°  nrocedure  given  in  Subroutine  FACTOR.  Subroutine  FACTOP  has 
been  retained  from  C0WCAC  I  fftT  factor  comnutatioric  not  ”tilizinrt  t^e  silhouette 
generator,. 

Subroutine  FACTOR 


This  subroutine  computes  configuration  factors  from  each  noint  on  Sur  ce  1 
selected  by  HAP  to  Surface  2.  The  exchange  coefficient  is  computed  by  numerical 
integration  of  configuration  factors  across  Surface  1,  from  which  the  form  fac¬ 
tor  is  finally  derived  as  the  area-weighted  mean  of  all  configuration  factors. 

Factors  are  computed  for  ea>.h  noint  along  each  manning  line,  movinp  from 
X-left  to  X-iight,  bv  translating  the  origin  of  the  Surface  2  coordinate  sys¬ 
tem  m  X.  The  analysis  and  equations  are  organized  for  minimum  conmutatxonal 
time;-  constants  at  each  loon  level  are  computed  once  prior  to  loon  entrv.  Be¬ 
cause  the  usual  outnut  desired  is  only  the  form  factor,  configuration  f-ctons 
ner  se  are  not  computed  unless  a  detailed  output  is  renuested.  \  numerical 
* itcgrstion  of  computed  noint  function  with  resoect  to  X  is  oeL*ormea  before 
proceeding  to  the  next  lino.  After  all  honzonal  integrations  are  rt  *nlcted, 
these  products  are  integrated  with  resn^ct  to  Y,  and  divided  bv  the  manning 
area  comnuted  m  subroutine  MAP., 
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A  standard  24  X  24  grid  results  m  62S  configuration  factors  to  be  :oro~ 
outed.  ’P'e  ouestion  naturally  arises  as  to  whether  this  nap"  configuration 
factors  are  actually  required.  If  the  configuration  factor  changes 
little  across  Surface  1,  then  it  is  nrobablv  too  many;  but  if  there  are  sham 
changes  in  the  factor,  and  third  place  accuracy  is  desired,  then  it  is  prob¬ 
ably  sufficient.  Contrary  to  off-hand  expectations,  a  more  sophisticated 
integration  rule  such  as  Simpson's  or  Weddle's  is  not  as  accurate  as  the 
trapezoidal  rule  for  standard  i"crc">e"rc  jf  the  factor  function  slope  changes 
raoidlv.  Weddle's  rule  was  initially  used  which  explains  why  the  program 
increment  control  is  m  groups  of  six  (except  the  initial  3  which  is  not  in 
ONFAC  I).  It  the  factor  varies  smoothly,  a  6  x  6  neddle's  rule  integration 
(49  factors)  is  probably  as  accurate  as  the  standard  62S  factors  nresently 
used  by  the  trapezoidal  rule.  The  time  saved  is  appreciable  when  running 
manv  factors.*  If  desired,  Weddle's  rule  may  be  inserted  in  the  source  deck 
and  compiled  with  no  other  changes  required. 

The  form  factor  computed  bv  the  above  is  from  that  part  of  Surface  1  which 
"sees"  Surface  2.  If  Surface  1  is  bisected,  then  the  computed  factor  must  be 
reduced  m  proportion  to  the  area  reduction.*  rhis  is  required  because  the 
total  active  side  of  Surface  1  entered  in  data  is  considered  the  radiant  sur¬ 
face., 

CENERAL  RULES  AND  RESTRICTIONS 

The  following  general  rules  and  restrictions  must  be  observed  for  normal 
progra<‘*  •  eration  • 

1.  Ail  '?”ta  must  be  derived  from  rignt-nandeci  rectangular  coordinate  svs- 

tens.,  “ 

2.  Points  1,  2  and  the  last  point  in  plane  surface  input  data  (Class  1 

and  4)  must  not  form  a  straight  line  in  space. 

3.  The  active  side  of  a  plane  or  nonnlanar  smface  is  established  bv  en¬ 

tering  the  boundary  points  in  counterclockwise  order,  as  thev  appear 
when  facing  the  active  side. 

4.  If  the  factor  to  a  Class  2  (nonnlanar)  surface  is  requuad,  onlv  the 
active  surfaces  should  be  seen  from  any  point  on  Surface  i ,  and  thev 
must  3 Iso  be  seen  from  every  point  on  Surface  1, 

5.  All  surfaces  used  as  Surface  2  which  utilizes  the  silhouette  generator 
(Classes  S,  6,  8,  or  4  if  included  m  Class  8)  must  npnerr  above  the 
plane  of  Surface  1,  i.e.,  all  Z  -  coordinates  must  be  noizero  and 
positive,  prior  to  factor  computation. 

6.  A  p'-iinar)  transformation  of  Class  8  data  is  not  permitted,  Also,  no 
auxiliary  transformation  is  permitted;  Surface  1  must  ,>e  in  tfe  XY 
plmk,  of  the  Multisurface  coordinate  system  as  entereu  in  data. 

7.  Detailed  restrictions  and  limitations  upon  ii  nit  data  are  given  m 
Section  IV. 
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srcrioN  iv 


INPUT  DATA 

DATA  SPECIFICATIONS  AND  SPECIFIC  RESTRICTIONS 

Innut  data  consists  of  externally  commuted  surface  data,  sneclfications 
for  intemallv  created  surfaces,  transformation  data  and  run  instructions 
vf'ctor  requests!.  Also  and  comments  cards  wav  he  ei.r^ed  as  required. 

Data  tync  is  classified  by  the  use  of  an  integer  fron  1  to  9  nl3ced  in 
coi'^n  2  of  the  data  nane  card,  followed  by  a  1  to  5  FORTRAN  character  name 
to  j>nmdc  data  identity  within  each  class.  The  classes  of  data  are  described 
below. 

ria<t<  l  -  Plane  Polygon 


The  X,  Y,  and  Z  coordinates  of  each  ooint  defining  the  surface  boundary 
arc  required.  Only  one  side  of  a  single  nlane  surface  can  be  nndo  active 
i,c,,  nnv  interchange  radiant  flux  with  another  surface.  The  active  side  is 
established  in  the  following  manner:  Face  or  look  at  the  desired  active  side, 
and  select  any  ooint  on  the  surface  boundarv  as  nomt  number  one.  Proceeding 
m  a  counterclockwise  direction  about  the  boundary  of  the  surface,  select  the 
remaining  noints  in  sequence.  If  this  rule  is  followed,  the  surface  will  al¬ 
ways  be  on  the  left  when  moving  along  the  boundary. 

The  Xf  Yf  and  Z  coordinates  of  each  ooint  arc  entered  on  the  data  cards 
in  the  above  seoucncc,  and  each  nomt  is  numbered  intemallv  according  to  its 
nosition  in  the  data. 

It  is  assumed  that  a  Class  1  surface  is  a  nlane  surface.  No  internal 
check  is  made  to  verify  this  (in  contrast  to  CONFAC  I).  If  a  substantial lv 
nonnlanar  surface  is  classed  as  a  nlane  surface,  serious  errors  in  manning 
could  result  if  used  as  Surface  1,  or  wrone  factors  cormutcd  if  used  as  Sur¬ 
face  2. 

No  ooint  connections  data  are  entered  under  Class  1;  the  silhouette  gener¬ 
ator  is  not  used. 

Class  2  -  Nonplanar  Surface 


Two  or  more  nlane  surfaces,  not  in  the  sane  nlane,  adjoining  or  connect¬ 
ed,  and  entered  as  one  nackagc  is  termed  a  nonnlanar  surface 

A  Class  2  surface  can  be  used  as  a  Surface  2  if  the  side  of  each  facet 
selected  as  the  active  side,  and  only  those  sides,  arc  seen  from  everywhere 
on  the  active  sido  of  Surface  1.  The  counterclockwise  order  of  data 
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entry  establishing  the  active  side  is  also  required  as  in  Class  1,  but  no 
orientation  vector  is  generated. 

No  connection  data  is  required  because  the  silhouette  peneratox  is  not 

used. 

Class  3  -  Internally  Generated  Plane  Polygon,  No  Connections  Data 

The  interna!  surface  generator  will  conmitc  the  coordinates  of  each  noint 
uc.'inmj*  a  nlanc  po'vyon,  oa:<sllcl  to  the  XV  olane,  with  ar.  orientation  vector 
erected  over  noint  1  and  directed  toward  the  +Z  axis.  A  detailed  description 
of  the  internal  surface  generator  is  given  m  Section  II. 

The  data  required  for  a  Class  3  surface  is: 

No.  of  sides,  3  -  N  1  100 

X  -  Axis  Radius 

Y  -  Axis  Radius 

X,  Y,  and  Z  coordinates  to  center  of  oolvpon 

A  Class  3  surface  is  used  in  the  same  wanner  as  a  Class  1  surface.  The 
sane  rules  and  restrictions  annly. 

Class  4  -  Plane  Polygon  with  Connections  Data  and  Class  5  -  Nonnlanar  Polygon 
or "Solid'  Surface  with  Connections  Data 


A  Class  4  surface  is  actually  a  Class  1  surface  with  connections  data 
added  waking  it  possible  for  it  to  be  processed  with  the  silhouette  renerator. 
Rut,  in  general,  no  useful  numosc  is  rained  hv  the  use  of  the 'silhouette  gen¬ 
erator  to  process  a  plane  surface,  unless  conbined  with  other  surfaces.  There¬ 
fore,  a  Class  4  surface  is  processed  as  a  Class  1  surface,  unless  it  is  listed 
under  a  Class  3  entry. 

A  Class  5  surface  is  always  processed  in  the  sinolc  node  l»v  SILFAC  unless 
listed  under  a  Class  S  entry. 

A  naxinun  of  100  boundary  points  wav  be  entered  describin'*  a  Class  4  or 
S  surface.  Uo  to  4  connecting  points  for  each  boundary  point  nr.v  be  entered. 

If  wore  than  four  connecting  points  are  required,  one  nav  enter  wore  boundarv 
data  Points  having  the  sane  coordinates  and  connecting  to  cadh  other,  using 
the  surplus  (3)  connections  to  satisfy  the  additional  connections  rcmiirencnt. 
However,  *  f  wore  than  two  such  identical  boundary  points  are  used,  the  surface 
cannot  be  processed  in  the  siwoYc  node”  This  restriction  in  rest  practical 
situations  can  be  circunvcntod  bv  separating  the  points  siifhtlv  with  little 
effect  on  the  final  forw  factor  comuted.  If  this  cannot  he  done,  the  surface 
rust  be  listed  under  r  Class  8  nanc  and  processed  in  the  cowpIcx  node. 


Class  6  -  Internally  Generated  Polygon  or  Polvhedyon,  Including  Connections 
Data  '  ' 


A  detailed  description  of  the  internal  surface  generator  is  given  m 
Section  II.  A  class  6  surface  is  always  processed  bv  SILFAC  -  in  the  simple 
node  if  used  directly,  and  conplex  if  listed  under  Class  8.  The  data  re- 
ouired  to  create  a  Class  6  surface  are: 

1.  No.  of  cress  sections 

2.  No,  of  cross  sections  divisions  (sides) 

3*  Coordinates  and  generating  wdn  of  first  cross  section, 

4.  2-coordinate  and  generating  radii  of  additional  cross  sections, 
if  any. 

Attention  is  directed  to  general  restrictions  4  and  5  in  Section  III, 

Class  7  -  Sphere 


The  radius  and  the  X,  Y,  and  Z  coordinates  of  the  sphere  are  renuired.  A 
primary  transformation  of  a  snhere  is  pointless  and  therefore  not  oemitted. 
Arbitrary  orientation  of  both  Surface  1  and  sphere  is  allowed.  One  peculiar¬ 
ity  exists  which  differs  fror  the  usual  treatment  of  bisected  surfaces.  If 
the  plane  of  Surface  1  bisects  the  sphere,  the  area  of  the  spherical  surface 
above  the  horizon  of  Surface  1  will  be  computed.  Now  the  sphere  cannot 
bisect  Surface  1  in  the  usual  sense,  but  it  is  Possible  a  bisected  snhere 
mav  be  partly  or  totally  inside  the  boundaries  of  Surface  1  -  embedded  in 
the  surface.  In  this  case,  the  program  will  merely  assign  a  zero  for  the  con¬ 
figuration  factor  when  the  viewpoint  from  Surface  1  is  inside  the  sphere.  This 
zero  will  ho  integrated  as  usual  with  the  other  factors  computed  alono  each  map¬ 
ping  line.  No  computation  of  the  Surface  1  area  seeing  the  sphere  is  made, 
however,  even  though  part  of  Surface  1  is  not  seen  bv  the  sphere. 

Class  8  -  Muitisurface 


A  Muitisurface  consists  of  from  one  to  eleven  Class  4,  5,  or  6  surfaces. 
A  Class  8  surface,  and  only  a  Class  8  surface,  is  processed  in  the  complex 
mode.  The  only  data  entry  necessary  to  indicate  the  surfaces  which  comprise 
a  Multisurfacc  are  the  names  assigned  each  surface. 

Class  9  -  Transformation  Data 


Transformation  Data  consists  of  the  coordinates  of  three  points  in  a  sur¬ 
face,  not  in  a  straight  line,  derived  t,om  the  ’’new  position  of  a  sin  face 
which  has  been  noved  in  its  coordinate  system.  One  nay,  with  coual  vnliditv, 
interpret  the  transformation  to  mean  that  the  origin  of  the  coordinate  svstem 
is  being  noved  to  a  different  position,  and  the  data  arc  the  coordinates  of 
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each  point  taken  fron  the  new  orioin.  The  three  oomts  selected  need  not 
he  chosen  or  entered  in  anv  particular  order,  nor  nust  the  sane  points  be  used 
if  norc  than  one  different  nr  mar''  transformation  oe  the  sane  surface  is  de¬ 
sire 

Run  Instructions 


Run  instructions  s«ecifv,  for  each  factor  desired,  the  foil  >nn<** 

1.  The  none  of  Surface  1  (cnittcr) 

2.  The  nare  of  Surface  2  (receiver) 

3.  Tronsfornarion  data  nane(s)  for  Surface  1  and/or  2,  i*  rcouxrcd, 

4.  Whether  a  standard  or  detailed  outnut  is  desired,  bv  inserting  code 
letter  "D"  for  detailed  outnut. 

5.  The  horizontal  and/or  vertical  divisions  to  he  used  m  nanpint*  sur¬ 
face  The  nator  divisions  which  nav  he  used  are  6,  12,  18,  24, 

30,  36,  42,  48,  54  and  60,  hut  m  run  instructions  these  divisions 
arc  specified,  respectively,  bv  the  integers  1,  2,  3,  4,  5,  6,  7,  8, 
9  and  10,  A  special  division  of  7  nav  he  snccificd  bv  the  integer 
11. 

DATA  OJMJ  >JSION  RISTRICTIONS 

1.  A  naxmun  of  100  houndarv  points  (300  coorunuttes)  for  each  surface 
entered  as  Class  1,  2,  4,  and  5. 

2.  A  naxmun  of  100  points,  enui valent  to  100  sides,  generated  bv 
Class  3  data. 

3.  For  Class  6  data,  the  nunber  of  sides  "lus  one,  tires  the  number 
of  cross  sections,  nust  not  exceed  101  if  nlnne,  and  102  if  non- 
nlanar, 

4.  The  f*rand  total  of  surfaces  entered  or  generated  bv  Class  4,  5 
and  6  must  not  exceed  11. 

5.  The  Rrnnd  total  of  surfaces  entered  or  rcncratcd  bv  Classes  1,  2, 

3,  4,  5  and  6  nust  not  exceed  26.  If  a  detailed  silhouette  output 
is  requested,  the  yrand  total  nust  not  exceed  16. 


6. 

The 

total 

nunber 

of 
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7  data  nust 
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not 
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PROGRAM  CONTROL 


The  program  deck  setup  is  shown  in  ripurc  36.  Note  the  Presence  of  a 
MT"  card  immediately  following  the  *  DATA  Card  and  the  variable  format.  A 
"T"  card  has  a  "T"  in  column  1,  and  serves  two  purposes.  Columns  2-72  nav 
contain  job  title,  name  of  programmer,  etc.  and  vail  be  nrinted  m  the  outnut 
of  innut  data  and  each  factor  result.  The  'T*  card  also  initializes  data 
Storage  locations,  so  that  new  innut  data  can  he  read  in.  This  means,  however, 
that  the  old  data  is  effectively  wined  out,  and  is  no  longer  available  for 
fac  or  commutations,  lulcs-  i«.*  entered  as  input  data.  *t  is  obviously  un¬ 
necessary  to  use  the  "T"  card  unless  all  available  locations  arc  used  un. 

Actiallv,  a  'T‘  card  docs  not  necessanlv  have  to  follow  the  variable  format 
unless  one  desires  the  title  to  be  nrinted,  because  the  data  location  counters 
arc  automatically  initialized  at  the  start  of  the  program,  But  subsequent  re¬ 
initialization  can  be  acconnlishcd  onlv  bv  a  card  with  a  "T"  in  column  1. 

It  appears  desirable  to  have  senarate  identification  of  the  various  fac¬ 
tors  connut''’1,  and  a  comments  card  has  been  provided  for  that  purpose,  The 
comments  card  has  a  "C”  in  column  1,  and  a  comment  nay  appear  in  columns  7  - 
72.  A  comments  card  nav  he  inserted  between  run  instruction  cards,  and  the 
uno  of  comment  given  on  the  card  will  be  printed  below  the  title  on  all  out¬ 
put  thereafter,  unless  superseded  b\  another  comments  card. 

Comments  outnut  nay  be  entirely  suppressed  bv  using  another  comments  card 
containing  blanks  in  columns  2  -  72. 

FORMAT 

All  data  may  be  entered  on  NAA  FORTRAN  Tixcd  10  Decimal  Data  sheets. 

Each  lino  represents  12  card  columns  with  six  lines  ner  card,  making  a  total 
of  72  card  columns  available  for  data  entry.  Columns  73  -  80  are  used  for 
card  identification  and/or  numerical  scoucr.cing  for  sorting  purposes. 

Title  Card 


A  title  card  is  characterized  by  an  alphabetical  "T"  placed  in  column  J. 
Columns  2-72  available  for  job  identification,  as  shown  on  Figure  9, 

Comments  Card 


A  comments  card  is  characterized  by  an  alphabetical  MC"  Placed  m  column 
1.  Columns  2-72  arc  available  for  run  comments,  as  shown  on  Firurc  f>. 

Surface  and  T -an s format ion  Data 


All  surface  and  transformation  data  is  preceded  by  a  name  card  uniquoly 
identifying  the  data.  A  nano  consists  of  six  FORTRAN  characters  (  a  connutcr 
"word”)  and  always  occupies  the  first  six  columns  of  the  nano  card.  The  data 
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class,  an  integer  fron  1  to  0,  oust  always  be  placed  m  colunn  1,  The  remin¬ 
der  of  t^c  nare  occupies  colunns  2-6,  and  it  is  important  to  note  that  a 
blank  smcc  is  considered  a  character  and  a  nart  of  the  nanc.  ror  cxanole, 
the  nanc  1S1  is  not  the  same  as  1  SI  or  1 _  SI, 

The  next  word  on  the  nanc  card,  colunns  7-12  nust  he  left  Man! .  The 
renainder  of  the  nanc  card,  (colunns  13  -  72)  wav  he  left  blank  or  core-cnts 
written  to  further  identify  the  data,  and  will  be  printed  out  alone  with  the 
nane  of  the  surface  as  nart  of  the  "Inout  Data"  nrint  out. 

The  dat i  identified  «»v  i„e  nane  card  oust  follow  the  n.mo  card.  Thcie 
arc  seven  different  fomats  which  nust  be  adhered  to  in  catering  data. 

Class  1  und  2 


The  nunber  of  noxnts  to  he  entered  describing  the  surface  appears  on  the 
first  line,  Figure  10,  followed  by  the  X,  Y,  and  Z  coordinates  of  each  noint 
in  sequence.  The  ord*r  in  which  the  points  arc  selected  in  the  surface  is 
explained  in  detail  in  Section  IV, 

Class  3 


The  nunber  of  sides  are  entered  on  the  first  line,  followed  bv  the  X,  Y, 
and  Z  coordinates  of  the  center  of  the  internnllv  qcncratcd  rolv<*on,  the  X- 
axis  radius  and  the  v-axis  radius  as  shown  »n  rimirc  11, 

Class  4 


The  total  nunber  of  points  describing  the  surface  are  entered  on  the 
first  line  as  shown  in  Figure  12.  The  X,  Y,  and  Z  coordinates  of  the  first 
noint  follow  on  the  next  three  lines,  The  fourth  line,  represent  me  12  col¬ 
unns,  is  divided  into  four  equal  parts  of  3  colunns  each,  Kach  noint  in  the 
surfaco  connecting  to  point  1  is  entered,  up  to  a  naxinun  of  four.  The  pattern 
is  repeated  for  the  renaininp  points  describing  the  surfaces. 

Class  6 


Tho  numbers  of  surface  cross  section  boundary  divisions  (sides)  is  r.ivcn 
on  tho  first  line  as  shown  in  Fipuro  13,  Tho  nunber  of  cross  sections  de¬ 
sired  is  specified  on  the  second  line,  followed  bv  the  X,  Y,  and  Z  coordinates 
of  the  base  (1st)  cross  section.  The  X-axis  radius  of  tho  base  cross  section 
is  pivon  on  the  last  lino  of  the  first  card.  The  Y-axis  radius  is  entered  on 
tho  first  line  of  the  second  card,  followed  bv,  if  j»oro  than  one  cross  section 
is  specified,  the  following,  repeated  for  each  crosi  section.  The  height  (Z- 
coordinate)  of  the  cross  section  above  the  XY  plane,  tho  X-axis  radius  and 
the  Y-axis, 
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Clas: 


7 


The  sphere  radium  i  entered  on  the  first  line,  followed  bv  the  X,  Y, 
and  Z  coordinates  locati***;  the  center  of  the  snhere  as  shown  in  Figure  14. 

Class  S 


The  nanes  of  the  surface  (s)  which  arc  to  be  entered  under  this  class 
arc  entered  together  on  one  card,  without  regard  to  order  as  shown  in  Figure 
15.  The  card  is  eouallv  divided  into  12  words  of  six  columns  each.  Fach  n<o>e 
.o  be  entered  mist  *» ;r  *  lenticsllv  in  n  word-snac*  as  it  appears  m  the 
word -space  on  the  data  name  card. 

Class  9 


The  first  noir.t  to  be  transformed  is  entered  on  the  rirst  line,  followed 
bv  the  X,  Y  and  2  coordinates  of  the  "new"  nosition  of  the  point  as  shown  in 
Figure  16.  The  second  nomt  to  be  transformed  immediately  follows  on  the 
fifth  line  followed  bv  th**  X-coordinatc  of  the  new  nosition  of  the  second 
point,  thus  completing  tn^  first  card.  The  Y  and  Z  coordinates  of  the  new 
position  of  the  sc  ond  point  arc  entered  on  the  first  two  lines  of  the  second 
card,  followed  by  the  number  of  the  third  point  to  be  transformed  and  its  new 
X,  Y,  and  Z  coordinates. 

All  of  the  numbers  entered  m  the  above  data  nnv  be  entered  as  fixed  or 
floating  point  numbers  except  connections  data,  which  must  be  entered  as  dec¬ 
imal  integers.  If  a  decimal  point  is  given  (fractional  nunbcis  must  have  dec¬ 
imal  points  viven),  the  floating  number  nav  be  located  anvwhcre  in  the  field 
(line);  if  no  decimal  point  is  given,  the  number  must  be  located  to  the  ex¬ 
treme  right  of  the  field  (no  blanks  to  the  right  of  the  number). 

Run  Instructions 


Six  FORTRAN  words  comprise  a  set  of  run  instructions;  two  sets  nay  he 
entered  on  one  card  as  shown  in  Figure  17.  The  first  set  starts  at  column  1 
and  the  second  set  starts  at  column  37.  Two  words  (12  columns)  comprise  one 
line  on  the  data  sheet.  The  name  of  the  Surface  1  data  is  entered  in  the 
first  word  (columns  1-6)  precisely  as  it  appears  in  the  first  word  of  the 
surface  data  name  card.  The  namo  o'*  the  Surface  2  data  is  entered  in  the 
second  word  (columns  3  -  12)  precisely  as  It  appears  in  the  first  word  of  the 
surface  data  name  card.  If  a  primary  transformation  of  Surface  1  is  desired, 
the  desired  transformation  data  name  is  entered  in  columns  13  -  18,  otherwise, 
it  is  left  blank.  If  a  primary  transformation  of  Surface  2  is  desired,  the 
name  of  the  transformation  data  is  entered  in  the  fourth  word,  ccltTO*  '9  - 
24.  If  a  detailed  output  is  desired,  the  alphabetical  character  "D"  is  entered 
in  column  2Zt  or  in  column  31.  If  the  "D"  appears  in  either  or  both  locations, 
a  detailed  output  will  result;  if  a  blank  is  in  both  locations,  a  standard  out¬ 
put  will  result.  Iho  horizontal  manning  division  "integer”  appears  in  column 
30,  unless  the  integer  10  or  11  is  used,  in  which  case  columns  29  and  30  arc 
utilized.  The  vertical  manning  division  "integer"  appears  in  column  36,  un- 
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FIGURE  14.  Clas3  7  Sphere  Specification  Input  Data  Format 
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FIGU  Z  16.  Class  9  Transformation  Data  input  Data  Format 


FORTRAN  FIXED  10  DIGIT  DECIMAL  DATA 


FIGURE  17.  Hun  Instructions  Input  Data  Format 
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less  the  integer  10  o*“  11  is  used,  in  which  case  coluans  3S  and  3o  arc 
utilized.  If  columns  29  and  30  are  both  blarl,  a  standard  integer  4  ncanine 
24  horizontal  divisions  of  each  nanmnp  line  is  used.  If  5S  and  36  arc  both 
blank,  a  standard  vertical  division  of  24  will  be  used*  The  above  fornnt  is 
rcncated  in  the  sane  nanner,  starting  from  colunn  37  on  the  fourth  line,  for 
the  second  set  of  run  instructions  on  the  card.  There  is  no  nuner i cal  linii 
to  the  nunher  of  run  instructions  which  nav  be  entered.  The  onlv  reouirenent 
is,  of  course,  that  the  data  called  for  has  been  loaded  in  under  the  nancs 
used. 
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SECTION  V 


PROGRAM  OUTPJrr 

Input  data  is  processed  ?pd  printed  out  prior  to  its  use  in  factor  com¬ 
putations  for  programmer  verification.  Toe  orientation  vector  head  end  is 
also  printed  out  for  all  plane  surfaces,  so  that  the  "active”  side  used  bv 
the  program  is  clearly  shown,  Class  3  and  6  specifications  as  read  in  are 
p> Sated,  along  with  the  surface  data  generated  r>v  the  specifications. 

A  standard  "aininun"  output  is  piven  when  the  code  letter  "D"  does  not 
appear  in  the  run  instructions,  consisting  of  the  following: 

1.  Run  number 

2.  Run  instructions 

3.  The  computed  form  factor 

4.  The  Surface  1  manning  area 

5.  The  exchange  coefficient  (/A  product) 

6.  The  total  area  of  Surface  1 

7.  If  Surface  l  is  bisected,  the  area  seen  hv  Surface  2 

8.  The  total  area  of  Surface  2,  if  Surface  2  area  can  be  computed 

9.  If  Surface  2  is  bisected,  the  area  seen  by  Surface  1,  if  that  area 

can  be  conouted 

10.  The  tine  in  seconds  snent  xn  Subroutine  SIEFAC,  if  utilised 

If  a  detailed  output  is  requested,  the  minimum  output  plus  the  following  is 
printed: 

1,  The  final  coordinates  of  Surface  1  and  Surface  2  prior  to  comnuta- 
tlon  of  configuration  factor*;, 

2,  The  X-Left  and  X-Right  coordinates  for  each  Y  division  of  Surface 
1  mapping,  including  horizontal  and  vertical  divisions  used, 

3,  Each  c on fi gyration  factor  computed.  Hie  output  is  g*v<*n  in  croups 
of  factors  easily  identified  because  the  last  factor  in  a  group 
occuoios, a  line  by  itself,  E.  ch  group  contains  the  configuration 
factors  computed  on  a  napping  line.  The  first  factor  in  the  group 
is  that  computed  at  Y-lcft  and  the  last  factor  in  the  group  is  that 


J  confuted  at  X-nght.  The  first  groun  renrcsents  the  first  manning 

5  line,  th»  second  groun  th^  second  manning  line,  etc. 

!  4.  If  the  silhouette  gcneiator  was  used,  the  silhouette  conoutcd  for 

'  points  selected  on  each  rapping  line  is  printed  out.  The  first 

mineral  given  is  the  manning  line,  the  second  is  the  ooint  on  the 
manning  line,  moving  fron  X-left  to  X-ripht.  The  numbers  following 
renresent  the  silhouette. 
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APPENDIX  A 


SAMPLE  PROBLEMS 


A  number  of  sa-ole  nroblens  have  been  devised  to  illustrate  the  capabil¬ 
ities  and  Imitations  of  CONPAC  II.  The  exannles  are  arranged  roughiv  in 
order  of  complexity,  bcemninp  with  simple  nlanc  surfaces  and  concluding 
wit!  a  complicated  "intervening  surface'*  nroblen  involving  nlanc  and  solid 
surfaces. 

The'  surface  configurations  unon  which  the  example  nroblens  are  based  are 
sh^wr.  m  acconoanvmc  illustrations.  Each  illustration  is  conveniently  group¬ 
ed  separately  with  the  nroblen  description  pertaining  to  the  surfaces  shown 
m  the  illustration,  along  with  the  input  data  sheets,  run  instructions,  pro¬ 
gram  outnut  and  a  short  discussion. 
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SAMPLE  PROBLEM  CROUP  A 


The  surfaces  shown  in  figure  18  are  similar  to  the  examples  riven  in 
CON FAC  I*  The  added  CON FAC  II  capability  of  bisectine  a  nonnlanar  surface 
is  demonstrated.  The  data  sheets  are  shown  m  Figure  19  and  the  results  are 
presented  in  Figure  20. 

Problem  l.\ 

In  Figure  18  (AI),  the  factor  between  t^e  rloor  of  x  x.u,  ical  room  (lFLOOPj 
and  an  adjacent  vail  (1*’AII.)  ;<•  confuted,  usin'*  standard  horizontal  and  verti¬ 
cal  nanomr  divisions  (24  x  24)  on  Surface  1.  A  detailed  output  is  rcoucstcd. 


Note  that  because  no  primary  ox  auxiliary  transformation  occurred,  the 
f '•''ordinate  system  is  the  same  as  the  innut  data  (unprioed)  coordinate 
system.  The  first  manning  line  starts  fit  the  origin  and  extends  to  noint  1 
in  1 FLOOR. 

Problem  2A 

In  Figure  18  (Al),  anv  plane  surface  may  be  used  as  Surface  1  nroviding 
it  has  been  nronerly  entered  in  data  prior  to  the  factor  renuest.  To  de¬ 
monstrate,  the  wall  (1WALL)  now  acts  as  Surface  1,  and  the  factor  to  the 
floor  (1 FLOOR)  is  requested. 

Note  that  Surface  1NAU.  is  not  in  tho  X-Y  plane  of  its  innut  (unnrimed) 
coordinate  system.  The  program,  therefore,  had  to  perform  an  auxiliary  transf¬ 
ormation  of  both  surfaces  to  the  priced  $’  stem  she~r.,  -rOr  tc  fac**'r  e ''Cita¬ 
tion,  to  get  Surface  1  in  the  XY  plane. 

Problem  3 A 

In  Figure  13  (A2j,  the  factor  from  the  floor  (1FL00R)  to  two  adiaccr.l  walls 
taken  together  (2HALLS)  is  requested.  This  is  a  valid  request  because  the 
boundary  data  describing  2 WALLS  form  a  valid  silhouette  of  2WALLS  from  anv 
noint  on  1  FLOOR.  The  factor  should  be  twice  that  to  one  wall  alone. 

Problem  4 A 

The  program  cannot  validly  compute  the  factor  from  a  nonnlanar  surface. 

A  Class  2  surface  is  assumed  nonnlanar.  The  factor  from  2NAILS  to  1 FLOOR  is 
requested  in  order  to  oiicit  the  diagnostic,  warning  the  user  of  mis  error. 

The  program  does  not  test  the  surface,  as  in  CONFAP  I.  If  a  nonnlanar 
surface  is  erroneously  entered  as  a  Class  1  surface,  it  will  not  he  rejected 
if  used  as  Surface  1  -  the  responsibility  lay  with  the  user  to  insure  that 
Surface  1  is  Planar, 
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Problcn  SA 


In  Figure  IS  (A3),  the  necessity  for  proper  order  in  data  entr/  is  empha¬ 
sized.  The  wall  data  are  deliberately  entered  in  a  clockwise  direction 
(1WALLR)  looking  at  the  active  surface,  instead  of  counterclockwise.  I,e"ee# 
the  orientation  vector  noints  in  the  wronv,  direction.  The  factor  fron  lFbOOR 
to  1WALLR  is  requested  in  order  to  elicit  the  diagnostic  which  alerts  the 
user  to  a  nossible  error. 

Prc1  len  6A 

In  Figure  is  (A*!),  CON'FAC  II  the  capability  of  bisecting  a  nonnlanar 
(Class  2)  surface.  The  factor  fron  irLOOR  to  2WALLZ ,  is  rcouested,  with  a 
tailed  output,  to  denonstrate  this  capability. 

Subroutine  bOICU  bisected  2”AILZ  at  the  XV  rlane,  and  reconstructed  the 
surface  by  eliminating  noints  2,  3,  and  4,  as  shown,  and  creating  new  noints 
2*f  3  ,  4  and  S  .  The  dashed  line  2*  3*  divides  Surface  1  (1TL00R)  into 
triangular  narts,  designated  A  and  B.  The  view  of  the  reconstructed  2KALI.Z 
fron  anvw^ne  in  a-ca  2  reflects  a  valid  silhouette  in  the  nroner  counter¬ 
clockwise  order.  When  reconstructed  2WALLZ  is  viewed  fron  area  A,  the  "o;nts 
still  fom  a  valid  silhouette,  but  the  order  is  reversed.  This  means  the 
computed  configuration  factor  will  be  to  the  hemispherical  snace  not  occupied 
by  2NAl!„Z,  and  will  be  negative.  So,  subroutine  FACTOR  subtracts  this  factor 
fron  1.0  to  yield  the  correct  factor  to  2NALLZ. 
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FIGURE  19.  Group  A  Sample  Problems  Input  Data  Code  Sheets 
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FIGURE  19.  Group  A  Sample  Problems  Input  Da^a  Cod©  Sheets 
(ecntlr.jicd) 
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FIGURE  19.  Group  A  Sample  Problems  Input  Data  Code  Sheets 
(continued) 
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I  FORTRAN  FIXED  tO  OGIT  OEOMAt.  DATA 


FIGURE  19.  Group  A  Sample  Problems  Input  Data  Code  Sheets 
(continued) 
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( continued) 


5 


(cor.tirued) 


000000000005000000CC('0000030000i,'''0900l}00000 
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( continued) 


FIGURE  20.  Group  A  Sanple  Problems  Program  Rcs’ilts 
(continued) 
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;roup  A  Sana,  Problo.-.s  Program  Kesui! 
(continued) 


(continued) 
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(contir.usc) 


FIGURE  20.  Group  A  Sample  Problems  Program  He suits 
(continued) 
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*ViPM.  ransim  r.«tnup  b 


* 

l 

i  The  fccnctrical  relationships  used  in  this  cxnnnle  arc  presented  in  ri<*- 

;  ure  21.  The  data  sheets  arc  shown  in  rieurc  22  with  results  in  Figure  23, 

Problcn  IB 

!  The  use  of  the  surface  generator  and  double  bisection  of  surfaces  is 

demonstrated.  The  "lane  Surface  1  PI. AT  I  is  entered  as  usual  in  the  data,  but 
»  the  octagonal  dirk  3PISK  js  cheated  bv  specifications  to  the  surface  genera¬ 

tor.  Koto  that  nc  connections  data  are  created  for  a  Class  A  surface,  but 
would  be  if  the  di«*k  were  nar.ed  CPISK. 

}  The  double  bisection  is  easily  seen  in  side  view  of  1PI.AT1  and  3DISK. 

|  The  results  of  tnc  factor  rcoucst  from  1PI.AT1  and  3DISK  is  shown  in  Run  #1 

outnut,  indicating  the  areas  in  each  surface  seen  hv  the  other.  The  number 
{  of  points  definin'*  3DISK  has  been  reduced  to  7  and  reorganized  because  of 

f  the  f>iscction,  as  seen  alonp.  the  dotted  line. 

| 

i  Problcn  2B 

Tlie  converse  factor,  3DISK  to  1PJ.ATJ,  is  requested  as  Run  #2,  Because 
the  dish  is  now  Surface  1,  the  final  coordinate  system  3PISK  is  aliened 

so  that  the  XY  nlaije^s  thc(nlanc  of  the  dish.  Point  1  becones  the  origin, 
and  line  segment  1  2  the  X  axis.  Note  that  the  exchange  coefficients  (/A) 
arc  verv  nearly  equal,  as  thev  should  be  because  of  the  reciprocity  theorem. 

Notice  that  the  factor  from  one  surface  to  the  other  nlonq  the  line  cf 
bisection  is,  in  realitv,  zero,  but  the  output  is,  in  some  cases,  non-zero 
though  quite  srall  (10’®  order  of  magnitude),  This  is  caused  hv  accumulated 
internal  truncation  error,  and  is  not  significant  enough  to  warrant  concern 
here,  (This  is  not  the  ease,  however,  with  so**c  silhouette  generator  con- 
nutations)  . 

\  Problcn  3B 

The  capability  of  coordinate  transformation  is  illustrated.  Run  #3  re¬ 
quests  the  factor  from  1PLAT1  to  3DISK  transformed  to  the  position  shown  hv 
the  transformation  data  9TDISK.  The  nropram  detected,  after  transforming 
30ISP,  that  it  bisected  1PLAT1,  As  the  cutout  shows,  tbo  part  of  1PI.AT1 
actually  manned  was  the  tranezoid  indicated  in  the  ton  View,  and  in  the  out¬ 
put  final  coordinate  data. 

Problem  4B 

it  Is  quite  feasible  to  jjenernto  or  manually  input  a  surface,  trai.sform 
tho  surface  to  a  different  location,  and  then  ask  for  the  factor  between  the 
original  surface  and  tho  transformed  surface.  This  is  shown  by  Run  *4, 
whero  3DISK  is  used  a*  Surface  1,  and  3PISK  transformed  by  9TDISK  is  used  as 
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Surface  2.  The  output  shows  a  bisection  of  3I)ISh,  removing  tie  4th  boundary 
point,  and  therefore  adding  a  roint  to  the  final  3IHSK  surface  boundaries, 
paking  it  9  instead  of  8, 

Problen  SB 

The  factor  fron  the  transformed  disk,  3DISK9TPICK,  to  1PLAT1  is  requested 
as  Kun  demonstrating  program  flexibility  in  that  Surface  1  is  now  trans- 
forned.  The  resulting  exchange  coefficient  is  very  ncarW  coral  to  Run  *5, 
as  it  should  be. 
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FIGURE  22.  Group  B  Sample  Problems  Input  Data  Code  Sheets 
(continued) 
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FIGURE  23.  Group  B  Sample  Problems  Program  Results 
(27  pages) 
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FIGURE  23.  Group  B  Sample  Problems  Program  Results 
(continued) 
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(continued) 
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FIGURE  2£.  Group  B  Sample  Problems  Program  :  e  >ult: 
(continued) 


102 


FIGURE  23.  Group  B  Staple  Problems  Program  Pesults 
(continued) 


OOOOOOOOOOOOOOOOOCOOOOO 


102 


FIGURE  23.  Group  E  Sar.ple  Problems  Progran  Results 
(continued) 
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FlGUliS  23.  Group  8  Sample  Problems  Program  Results 
(continued) 


FIGURE  2.3.  Group  B  Sample  Problems  Program  Results 
(continued) 


FIGURE  23.  Group  B  Sanpie  Problems  Progran  Result? 
(continued) 
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FIGURE  2 3*  Group  B  Sample  Problems  Program  Results 
(continued) 
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(continued) 
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FIGURE  23.  Group  B  Sample  Problems  Program  Results 
(continue) 


Ill 


conti 


FIGURE  23.  Group  B  Sample  Problems  Program  Results 
(continued) 
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FIGUH2  23*  Group  3  Sacple  Problem  Program  Results 
(continued) 


FIGURE  23.  Group  B  Sample  Problems  Program  Re: 

(continued) 


FIGURE  23*  Group  B  Sample  Problems  Program  Results 
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FIGURE  23-  Group  B  Sample  Problems  Program  Results 
(continued) 


jfcJI  ^<^l*s**v**, ******  ■*' 


s 


I 

i 


SAMPLE  PROBLEM  GROUP  C 

The  geometrical  relationships  for  this  sannlc  nroblen  are  shown  in  Fieure 
24,  The  data  sheets  are  presented  in  Figure  25  and  the  results  arc  shown  in 
Figure  26. 

Problen  1C 

In  this  nroblc**,  a  solid  surface  which  could  not  be  created  bv  the  nroeran 
su?*«  e  generator  is  <J».Ureu  H&niKiUy  along  t,itli  the  necessary  connections 
d-tn.  A  cube  with  *our  truncated  comers,  na red  SO*BF,  s<  entered  in  data 
fron  a  convenient  location  in  its  coordinate  svsten,  i.e.,  at  the  origin,  as 
shown  in  Figure  24.  Only  three  points  were  computed  and  entered  as  9TCUBh 
transformation  data  to  move  the  surface  to  the  desired  position  shown  over 
Surface  1,  1PLAT5.  The  factor  fron  1PLATS  to  5CUBE9TCHBE  is  requested  as  Run 
#1. 


The  silhouette  generator  was  used  to  compute  the  silhouette  rrom  each 
point  in  1PLAT5,  and  because  a  detailed  output  was  requested  with  G  horizontal 
and  6  vertical  divisions  of  1PL\TS,  49  silhouettes  were  computed  as  shorn  in 
Figure  26.  The  nmbers  following  each  identifvinr  naooine  line  and  manning 
point  number  arc  the  boundary  point  numbers  which  form  the  silhouette  when 
connected  together.  It  was  possible,  since  there  are  no  crossovers  in  the 
silhouette,  to  run  this  nroblcn  in  the  simple  node  in  SIIF/C  at  "renter  speed. 
The  rapping  divisions  were  deliberately  set  at  6  x  6  to  reduce  the  cutout, 
some  exncrinentation  is  renuired  to  detemine  how  nnny  divisions  are  renuired 
to  yield  the  factor  to  the  accuracy  desired. 

Problem  2C 

fhe  silhouette  generator  requires  that  all  points  in  Surface  2  he  above 
the  plane  of  Surface  l  when  operating  in  either  the  simple  or  complex  mode, 

A  view  of  5CUBK  in  its  original  position  fron  1PLAT5  clearly  shows  part  of 
SCUBK  below  the  surface  of  1PLAT5;  the  run  is  therefore  rejected  with  a  diag¬ 
nostic  indicating  this  condition. 
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FIGURE  25.-  Group  C  Sample  Probleos  Input  Data  Codn  Sheets 
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FIGURE  26.  Group  C  Sample  Problems  Program  Re suit ?* 
(7  pages) 


FIGURE  26.  Group  C  Sample  Problems  Program  Results- 
(continued) 
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FIGUFwB  26.  Group  C  Sample  Problems  Program  Results 
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FIGURE  26.  Group  C  Sample  Problems  Program  Results 
(continued) 
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The  geometrical  relationship  for  this  sanole  problem  are  presented  in  Fig¬ 
ure  27.  The  data  sheets  are  shown  in  Figure  2S  and  the  results  are  presented 
in  Figure  29. 

Problen  ID 


The  rcfcrrrccd  figure  27,  shows  a  truncnted-caro-on-cylinder  and 

a  disk,  skewed  with  resncct  to  the  cylinder-cone  centerline.  The  cylinder- 
cone  is  created  by  the  surface  generator  as  6CYTR,  a  32-side  1  solid  in  its 
final  position  in  the  unprined  coordinate  system.  The  disk  is  also  internally 
generated,  but  because  the  generator  (in  its  nresent  version)  is  limited  to 
cross  sections  parallel  to  XY  nlane,  the  disk,  SilSKC,  had  to  he  transformed 
to  the  skewed  position  by  transformation  data  9TDSKC.  The  results  arc  shown 
in  Figure  29.  The  simple  mode  was  used  for  processing  because  no  line  seg¬ 
ment  crossovers  arc  nresent,  which  enabled  the  use  of  transformations  to 
construct  the  nroblem.  The  warning  note  concerning  the  difference  between 
the  napping  area  and  the  actual  3DISK  surface  area  is  suppI led  to  attract  at¬ 
tention  to  possible  errors  in  Surface  1  data  entry  of  the  choice  of  mapping 
increments.  As  indicated  in  the  comments  on  Problen  1  of  Sample  Problem 
Group  C,  the  coarse  increment  6x6  was  selected  to  reduce  output.  A  finer 
increment  should  probably  be  used  to  insure  accuracy  to  the  third  place,  if 
such  is  desired.  It  must  be  emphasized  that  the  form  factor  obtained  in  Run 
*1  is  tho  factor  to  the  solid  figure,  6CYLTR,  which,  of  course,  includes  the 
bases.  Since  the  factor  to  the  skin  is  the  desired  number,  it  is  ncccssarv 
to  subtract  the  factors  to  the  ends.  The  upper  end  is  obviously  not  seen 
(/  =0),  so  the  factor  to  the  base  only  must  be  obtained.  The  base  is  easilv 
created  by  the  surface  generator  (3DISKB),  but  it  is  created  with  the  orienta¬ 
tion  vector  pointing  toward  tbe  *Z  axis--the  wrong  way.  It  is  necessary  to 
turn  it  around  by  a  primary  transformation-- (9TDSKD).  Thus,  the  full  capabil¬ 
ity  of  the  nrinary  transformation  feature  js  utilized  and  exemplified,  show 
by  Run  #2.  The  factor  to  the  skin  of  6CYLTR  is  obtained  b”  subtracting  the 
results  of  Run  *2  from  Run  #1,  or 

■^skin  "  ^total  "  ^basc 

7  -  0.18946  -  0,09955 

/  «  .08991 

The  exchange  coefficient  is  computed  in  a  similar  manner. 
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FIGURE  28.  Group  D  Sample  Problems  Input  Data  Code  Sheets 
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FIGURE  29.  Group  D  Sample  Problems  Program  Results 
(14  pages) 


FIGURE  29.  Group  D  Sample  Problems  Program  Results 
(continued) 
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FIGURE  2$.  Group  D  Sample  Problems  Program  Results 
(continued) 
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(continued) 
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FIGURE  29-  Group  D  o&mple  Problems  Program  Results 
(continued) 
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(continued) 
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FIGURE  29.  Group  D  Sample  Problems  Program  Results 
(continued) 
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SAMPLC  PROBLEM  CROUP  E 


The  capability  of  obtaining  form  factors  to  spheres  in  any  rvositicn  rel¬ 
ative  to  Surface  1  is  demonstrated,  Figure  30.  Closed  form  configuration 
factor  solutions  aro  utilized,  enabling  very  ranid  comnu tat ions.  The  data 
sheets  are  presented  in  Figure  31  and  the  results  in  Figure  32. 

Problem  IE 

The  factor  from  a  rectangle,  1PLAT8,  to  a  snhere,  7SPII),  fullv  above  the 
plane  of  1PLAT8  (Case  I)  is  Yeoucsted  as  Run  #1.  The  configuration  factor 
solution  in  this  case  is  extremely  sinole  (see  Apnendix  E),  which,  in  addition 
to  the  coarse  canning  of  1PLAT8,  accounts  for  the  short  computational  time. 

Problem  2E 

The  factor  from  1PLAT8  to  7SPH2  is  reoucsted.  The  snhere  is  the  sam,  size 
as  7SP1I1,  except  part  or  the  lower  half  of  the  snhere  is  below,  and  outside 
of,  the  surface  of  1PLAT8  (Case  II).  The  results  are  shown  as  Run  #2. 

Problem  3E 

A  larger  sphere,  7SPI13,  is  located  with  mart  of  the  unner  half  of  tho 
sphere  below  the  surface  of  1PLAT8  (Case  III),  The  results  are  shown  as  Run 

#3. 

Problem  4E 

The  nrogram  will  also  commute  the  factor  to  a  snhere  which  is  embedded 
in  Surface  1,  illustrated  bv  7SPII2  and  1PLAT7.  However,  no  attennt  is  nade 
to  determine  what  nart  of  1PLAT7  sees  tho  snhere;  when  «  manning  noint  on 
Surface  1  anpears  inside  tho  sphere,  a  configuration  factor  of  zero  is  re¬ 
turned  and  integrated  along  with  the  other  commuted  factors.  Therefore,  in 
Run  *4  we  s^e  no  indication  that  1PLAT7  is  bisected  by  7SPH2,  although  in 
reality  it  is.  The  nroblem  is  handled  in  this  wav  because  of  ttfe  extreme 
complexity  of  the  general  determination  of  thai  mart  of  Surface  1  not  seen 
by  the  sphere. 

Problem  SF, 

Tho  trivial  caso  of  the  snhere  completely  below  Surface  i  is  illustrated 
by  Run  *5. 
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FIGURE  30,  Group  E  Sample  ProbKvts  Geometry 
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FIGvPJS  31.  Group  E  Sanple  Prohleas  Input  Data  Codo  Sheets 
(continued) 
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FIGURE  31.  Group  E  SanpTe  Probiers  Input  Data  Code  Sheets 
(continued) 
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FIGURE  32.  Group  S  Sample  Problcr.3  Program  Results 
(13  pages) 


(continued) 


FIGURE  32.  Group  E  Sample  Problems  Program  Results 
(continued) 


FIGURE  32.  Group  2  Sample  Problems  Program  Results 
(continued) 
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COORDINATES  OF  POINTS  ON  »OONOA*Y  OF  SURF 
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FIGURE  32.  Group  2  Sample  Problems  Pi-ogran  Recult;; 
(continued) 
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(continued) 
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(continued) 
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F2GUB3  32.  Group  S  Sample  Problems  Program  Results 
(continued) 


rlGURS  32.  Group  E  Sample  Problems  Program  Results 
(continued) 


Figure  32.  Group  E  Sample  Problems  Program  Results 
(continued) 
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Figure  32.  Group  E  Sample  Problems  Program  Results 
(continued) 


SAMPLE  PROBLEM  GROUP  F 


The  capability  of  computing  factors  to  surfaces  which  are  occluded  bv 
intervening  surfaces  is  demonstrated  as  shown  in  Figure  33.  The  data  sheets 
are  presented  in  Figure  34  and  the  results  in  Figure  35. 

Problem  IF 

The  factor  from  1PLA10  to  1PLA9  is  desired  with  surface  1PLAU  interven¬ 
ing.  The  surface  SCOPLA.  representing  the  boundary  points  of  1PLA9  including 
corrections,  and  the  boundary  points  of  lPLAii  including  co^cttions  is  enter¬ 
ed  in  data.  Because  the  silhouette  is  complex,  the  surface  SCOPLA  must  be 
reentered  as  SCOPLA  to  enable  the  silhouette  generator  to  operate  in  the  com- 
nlex  mode.  The  factor  to  8C0PLA  from  jPLAIO  is  requrested  as  Run  #1.  Then., 
the  factor  from  1FLA10  to  1PLA11  is  requested  as  Run  *2,  The  factor  from 
1PLA10  to  1PLA9  is  merely  the  difference  between  the  two, 

/-  0.26787  -  0.20146 

f  m  0.06641 


Problem  2F 

This  problem  also  illustrates  the  capability  of  determining  factors  to 
occluded  surfaces,  but  data  is  entered  and  handled  in  a  different  manner. 

The  factor  from  1PLA10  to  6PIPE2  is  desired,  taking  into  account  the  flux 
intercepted  by  SPIPE1. 

The  coordinates  defining  6PIPE2  are  internally  generated.  5PIPE1  is  en¬ 
tered  manually,  and  the  two  surfaces  are  combined  for  complex  processing  as 
820FCM. 

Notice  that  SPIPE1  includes  a  line  segment--a  ’’bridge"  line— connecting 
point  7  in  5P1PE1  to  6PIPE2.  If  this  line  or  any  other  suitably  oriented 
line  serving  the  purpose  were  not  present,  then  the  silhouette  generator 
would  not  include  6PIPE2  in  any  of  the  silhouettes  computed  from  points  on 
mapping  lines  6  and  7  on  1PLAI0.  The  line  does  not  have  'o  actually  be  in 
any  surface— it  need  only  appear  to  intersect  both  surfaces  in  the  silhouette. 

The  form  factor  to  820FD1  is  0.28139  (Run  #3);  to  5PIPE1  alone  is  0,21556 
(Run  #4);  therefore,  the  form  factor  to  6PIPE2  is  the  difference  or  0,06583, 

Problem  3F 

This  problem  Illustrates  improper  use  of  the  program,  and  ir  Particular, 
the  silhouette  generator.  The  factor  form  1PLA12  to  SCOPLA  1*  request'd  as 
Run  *5.  Note  that  the  data  SCOPLA  is  In  eallty  two  surfaces.  These  surfaces 
when  viewed  from  1PLA10  or  1PLA12  present  a  complex  silhouette,  and  therefore, 
must  be  processed  In  the  complex  mode.  However,  when  a  class  4,  5  or  6 
surface  is  specified  as  Surface  2,  the  simple  mode  is  always  used.  The 
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silhouette  generator  conseouer.tly  saw  only  1PLA11  sometimes  and  only  1PLA9 
sometimes;  this  condition  would  not  be  relieved  bv  use  o*  a  bridge  line,  be¬ 
cause  the  total  silhouette  is  complex,  uqd  *vust  be  nrocessed  as  complex. 

Problem  4F 

When  a  class  8  surface  is  used  Surface  2,  Surface  1  must  be  in  the  XY 
plane  of  the  Surface  2  coordinate  System,  with  its  orientation  vector  pointing 
toward  the  *Z  axis.  The  rcsutt$;£?  Run  #6  show  the  diagnostic  resulting  from 
a  request  for  the  factor  from  1PLA12  to  3C0PLA. 
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Figure  34*  Group  F  Sample  Prcblosa  Input  Date  Code  Sheets 
(continued) 
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Figure  34.  Croup  F  Sample  Problems  Input  Data  Code  Sheets 
(continued) 
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Figure  34.  Group  F  Sanpla  Problems  Input  Data  Code  Sheets 
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Figure  34.  Group  F  Cample  Problems  Input  Data  Code  Sheets 
(continued) 
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Figuro  35,  Group  F  Sample  Problems  Program  Results 
(18  pages) 


Figure  35.  Group  F  Staple  Problems  Program  Ke suits 
(continued) 


Figure  35.  Group  F  Sample  Problems  Program  Results 
(continued) 
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Figure  35.  Group  F  Sample  Problems  Program  Ke&ult s 
(continued) 
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Figure  35-  Group  F  Sample  Problems  Program  Results 
(continued) 


Figure  35.  Group  F  o ample  Fro  bio.  as  Program  Results 


Figure  35.  Group  F  Sample  Problems  Program  Rooulto 
(continued) 


figure  35*  Group  P  Sample  Problems  Program  Results 
(continued) 


■55’’' 


}' 

fe 

I* 


§ 

t 


r 

f 

* 

* 

i 


i 

\ 


183 


’rogr.'tn  Results 


Figure  35.  Croup  F  Sample  Problems  Program  Results 
(continued) 


186 


AMTXDIX  B 


PROGRAM  DrCK  srrup,  listings,  and  maps 


The  nropran  deck  arrangement  shown  in  Figure  36  contains  a  nain  nrogram 
and  six  subnroprams  which  are  listed  in  this  annendix.  A  listing  of  the  nain 
nrogran,  7J360,  is  shown  in  Figure  37  followed  by  a  nan  of  the  core  storage 
1  ications  in  Figure  3S« 

The  first  subprogram  in  the  deck  setun  is  subroutine  IWIVT.C  which  is 
shown  in  ripurc  39  and  the  nan  of  core  storage  in  Figure  40.  The  transforna- 
tion  subroutine,  TXCRM,  is  nresented  in  Figure  41  and  the  nan  of  core  storage 
in  rigure  42.  The  listing  .and  man  of  subroutine  DOICti  is  nresented  in  rigurcs 
43  and  44.  The  listing  and  nan  of  subroutine  ft\P,  is  nresented  in  Figure  4S 
and  46.  Thc  listing  of  the  subroutine  FACTOR  is  nresented  in  Fipurc  47  and 
the  nan  of  *  "  -^nrace  ir  Fieure  4S.  Subroutine  SiLFAG  listing  and  core 
storage  is  1  ure  SI  shows  the  variable  fornats 

used  bv  thi 

This  IBM  FORTRAN  II  nrogram  uses  two  innut-outnut  statenents  which  nust 
be  modified  for  computing  svstens  other  than  the  NAA  flonitor,  FIB  III,  svsten. 
These  arc 

READ  N,  List 
PRINT  N,  List 

A  convenient  FAP  assembled  prop. ran  is  included  which  will  convert  the 
READ- PR I NT  statement  to. 


READ  INPUT  TAPC  S,  N,  List 
WRITE  OlfTPUT  TAPE  6,  N,  List 

This  nsscnbly  is  listed  m  Figure  52.  The  convert  to  any  other  commuting 
svsten  using  norinhcral  enuionent  and  not  using  the  sane  tane  designations, 
the  last  three  toil  cards  arc  sinnlv  changed  to  read 

MIN  LOU  A  (Innut  statement  tame  number) 

MOlfT  EOU  B  (Output  statenent  tane  number) 

HPUNOf  EOU  C  (runch  statement  tane  nurfcir) 

Tor  conouting  centers  using  attached  printing  enuinnent,  the  PAP  asscnblv 
can  be  renoved  and  the  program  will  execute  in  that  system. 
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APPENDIX  C 


COORDINATE  TRANSFORMATION 


PRIMARY  TRANSFORMATION 

As  indicated  in  Section  II,  the  surface  coordinate  transformation 
technique  employed  by  the  program  does  not  require  transformation  para¬ 
meters  such  as  direction  cosines,  Euler  angles  and  translation  terms  to  be 
o.'terea  as  inpu;  aa-a  i-nc  orogr&za  ror  transiormation  purposes.  Instead, 
the  x,  y  ani  z  coordinates  of  three  points,  not  in  a  line,  ari  given  from 
the  new  origin .  or  to  the  new  position  of  the  surface.  These  data  are 
then  used  to  derive  the  rotational  and  translation  terms  required  to  trans^ 
fora  tne  remaining  surface  data  to  the  new  origin  or  surface  position. 

The  classical  equations  lor  transformation  of  rectangular  coordinates 
in  space  are  employed  for  both  primary  and  auxiliary  transforma i ion.  The 
new  x,  y  and  s  coordinates  in  terms  of  the  old  coordinates  are: 

x  ■  x'cos  +  y’cos  C’^  +  z'cos  Qtj  +  R  (1) 

y  *  x’cos  &i  +  y'cos  82  +  z*cos  +  L  (2) 

z  ■  x'cos  Y^  +  y'cos  y2  +  z’cos  y^  +  K  (3) 

Note  that  there  are  9  unknown  direction  cosines  and  3  translation  terms,  or 

a  total  of  12  unknowns.  It  is  clear  that  the  coordinates  of  four  points 
from  the  new  origin  are  required  if  these  equations  are  to  be  used  directly 
to  determine  the  unknown  parameters. 

It  can  be  shown,  however,  that  the  coordinates  of  three  points  (not 
in  a  line)  are  sufficient  and  necessary  to  fix  the  position  of  a  surface 
in  any  rectangular  coordinate  system.  It  appears,  therefore,  that  another 
point  must  be  made  available  for  solution  of  the  above  equations,  or 
another  technique  developed  which  directly  requires  only  three  points. 
Investigation  of  the  latter  yielded  a  complex,  difficult  to  program,  solu¬ 
tion.  On  the  other  hand,  solution  of  the  classical  equations  is  straight¬ 
forward,  but  requires  the  extra  point  (not  in  the  plane  of  the  other 
threo).  Rather  than  require  the  user  to  supply  the  extra  point  in  data, 
it  was  decided  to  generate  the  point  as  a  unit  normal  vector  above  the 
second  point  given  in  transformation  data.  This  extra  point  must,  of 
course,  bo  generated  in  both  old  ani  now  coordinate  systems. 

Figure  53  depicts  a  primary  transformation  or  Surface  A  from  the  old 
(primed;  to  tnc  new  (unpriced)  coordinate  system.  Note  that  the  pr^ary 
transformation  shown  affects  only  one  urface,  whereas  both  surfaces  ar<? 
involved  ainr'iMury  will  be  disCUSSOd  in  ©Ore 
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Figure  53.  Coordinate  Transfomation 


><  * 


detail  later  in  this  Appendix.  The  angles  (o',  0  and  y)  shown  in  Equations 
1,  2  and  3  are  related  in  the  following  manner  with  the  primed  and  unprimed 
coordinate  axes  shown  in  Figure  53: 

«  ZO'X'OX,  a2  -  A)'Y'a,  or3  «  /0*2'0X 

%  «/o'X‘or,  e2  *=■  dO’Y'ci,  e3  «zcz •or 

Y,  “^b'X'OZ,  v2  "d’OTOZ,  y3  »/0>Z*02 

Given  the  coordinates  of  points  2,  5  and  6  in  Surface,  A,  and  tne  gen¬ 
erated  point  U,  iron  both  the  old  and  new  coordinate  system/.,  we  may  write 
four  independent  equations  similar  to  Equation  1#  Using  Equation  1,  the 


resulting  set  of  equations  in  x  is : 

x>  •  xj  cos  «1  +  y2  co«  a2  +  zi  cos  +  11  (4) 

x<;  *  x^  cos  Of,  +  y^  cos  »2  +  z£  cos  03  +  H  (5) 

x6  *  cos  °1  +  cos  &2  +  zl  cos  °3  +  H  (6) 

Xy  -  x^  cos  otj  +  yr*  cos  <xz  +  z ^  cos  o^  +  H  (7) 


We  may  similarly  write  two  more  sets  of  equations  in  y  and  z  similar  to 
Equations  2  and  3,  for  a  total  of  twelve  independent  equations*  Each  set 
of  four  simultaneous  equations  is  solved  by  Cramer's  Rule  (Reference  2) 
for  the  unknown  direction  cosines  relating  the  old  y*  and  z'  axes  to  the 
new  x,  y  and  z  axes. 


For  example,  using  the  set  developed  above  for  Equation  1  (Equations 
A,  5,  6  and  7),  the  coefficient  determinant  D  is 


xl 

*u 


75 

n 


1 

1 

1 

1 


(8) 


By  Cramer's  Rule,  we  successively  replace  the  elements  in  each  cdlxian  of 
the  set  with  the  respective  element  on  the  left  of  each  equation  in  the 
set.  For  example,  the  solution  for  coc  ct  is: 


271 


5 

x6 

*U 


(9) 


a.xi  similarly, 


y2 

y5 

y6 

yu 


*2 

x5 

x6 

*U 


(10) 


The  above  process  is  repeated  for  the  solution  of  cos  EL,  cos 
cos  Y2*  ejr-d.  cos  Y3#  using  the  sets  developed  for  y  and  z.  Tie  coefficient 
determinant  is  the  s<w*j  * Or  all.  sets,  because  the  coefficients  of  the  un¬ 
knowns  in  all  sets  are  identical, 

To  increase  computational  efficiency,  repeated  factors  in  the  expanded 
determinants  are  computed  only  once  for  each  set.  Also,  considerable  econ- 
ony  results  in  computing  the  unknown  direction  cosines  cos  oj,  cos  cos  XL, 
as  the  cross  product  of  the  corresponding  direction  cosines  of  the  other 
(y  and  z)  axes: 


COS  0^  "  COS  @2  C03  Y3  “  COS  Y2  COS 

(n) 

cos  *■  cos  Y2  cos  03  -  cos  <*2  005  Y^ 

(12) 

cos  Yj,  “  cos  a2  c°8  ^3  ~  008  ^2  c°8  al 

(13) 

The  translation  components  H,  L  and  K  are  computed  by  substituting  the  co¬ 
ordinates  »of  the  point  below  point  U  in  the  surface  to  be  transformed  into 
Equations  1,  2  and  3,  along  with  the  known  values  of  direction  cosines. 

For  the  Surface  A  shewn  in  Figure  53,  using  point  5, 


Xj  cos  <*2  -  yj  cos  ■  z5  005  ^3 

(1A) 

x5  C08  cos  B2  “  c05  “3 

(15) 

x5  008  Y1  "  y5  c05  Xs  ”  *5  cos  T3 

(16) 
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The  program  now  transforms  all  point  coordinates  in  Surface  A  from  the  old 
to  the  new  system  by  direct  substitution  in  Equations  1,  2  and  3. 

The  method  outlined  above  will  always  perform  the  transformation  de¬ 
sired,  providing  the  three  points  selected  are  (1)  sufficiently  separated 
in  space,  (2)  accurately  computed,  and  (j)  do  not  form  a  straight  line  in 
space.  Because  the  fourth  point  U  is  always  computed  outside  the  plane  of 
the  other  three,  the  coefficient  determinant  D  can  never  be  0.  Hence,  by 
Cramer’s  Rule,  a  unique  solution  must  always  exist. 


AUXILIARY  TRAN SFOPJKAT ION 

An  auxiliary  transformation  transforms  the  coordinates  of  both  sur¬ 
faces  into  the  reference  plane  of  a  specified  surface — the  "Control"  sir- 
face.  In  Figure  53,  Surface  A  is  the  control  surface;  the  auxiliary  trans¬ 
formation  depicted  transforms  both  Surface  A  and  B  from  the  unprimed  (old) 
system  to  the  double-primed  (new)  system.  In  general,  the  origin  0"  in 
the  control  surface  Is  always  point  1  in  the  control  surface  coordinate 
array.  It  may  not  always  be  the  first  point  entered  in  input  surface  data; 
if  a  bisection  of  the  surface  occurs,  and  the  original  point  1  is  not  seen 
W  Ihe  ulher  3.uTac<j  (assumed  planar),  then  a  new  point  1  will  be  computed. 
The  new  point  will  be  used  as  the  origin  0".  The  same  processing  occurs 
for  internally  generated  surfaces.  Only  surfaces  classed  as  plane  sur¬ 
faces  may  bo  control  surfaces,  i.e.,  Classes  1,  3,  4  and  planar  6.  For 
example,  if  both  Surface  A  and  B  in  Figure  53  are  plane  and  bisect  each 
other,  t>£0  auxiliary  transformations  would  occur  to  facilitate  surface  re¬ 
construction.  Actually,  if  Surface  A  were  entered  as  Surface  1,  the  first 
auxiliary  transformation  to  occur  would  be  to  point  1  in  Surface  B,  rather 
than  Surface  A,  Ttiis  would  not  occur,  however,  if  Surface  A  were  po£.  bi¬ 
sected  by  Surface  B.  In  any  case,  the  last  transformation  always  is  to 
point  1  in  Surface  1,  so  that  mapping  and  factor  computation  may  proceed 
forthwith. 

The  processing  of  an  auxiliary  transformation  differs  from  the  primary 
transformation  because  unknowns  may  be  more  readily  computed  from  available 
data.  Equations  1,  2  and  3  are  rewritten  for  the  auxiliary  old  and  new  co¬ 
ordinate  systems. 


x"  ■  x  cos  +  y  ecs  z  cos  +  H"  (17) 

yM  -  x  cos  8^  +  y  cos  S2  +  z  c05  83  +  L"  (18) 

z"  “  x  cos  +  y  cos  y 2  +  z  cos  Y^  +  K"  (19) 

The' angles  ure  defined  as  follows,  referring  to  Figure  53: 

-ZQXOVX",  a2  YOYO’S",  Oj  «/0Z0"X" 
bj  -ZOXO"!*,  B2  “Z0Y0"Yn,  &3  -Z0Z0"Y" 

Yj,  -  ^0X0*2",  Y2  -^0Y0"Z",  Y3  -^0Z0"Z" 
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Because  the  0"XJ:  axi3  in  the  new  system  is  directed  along  the  line  segment 
formed  by  the  first  and  second  point  ir.  the  control  surface,  the  direction 
cosines  related  to  that  axis  are  readily  computed. 

The  length  of  line  12  in  Surface  A  is 

1&12  ^*2  "  +  *y2  “  yl'2  +  ^zZ  ~  zl?  > 

and  the  direction  cosines  relating  the  new  0"XM  axis  to  the  old  OX,  OX  and 
OZ  axes  are: 


cos  o^  “  (x2  -  x1)/LS12 

(20) 

ccs  a2  “  (y2  -  yi 

(21) 

cos  *3  =  <2o  -  ^y)/^i2 

(22) 

Because  th»,  new  OnZ‘'  axis  is  directed  alc:g  tne  surface  unit  or: 
tor  (point  0  above  Surface  A)  the  cosines  relating  that  axis  to 
OX,  OY  and  OZ  axes  are 

Lcntaticn  vec- 
the  old 

008  Yl  -  *b  “  *1 

(23) 

cos  y2  -  y0  -  7l 

(24) 

cos  V3  “  zQ  “  21 

(25) 

The  remaining  direction  cosines  are  computed  by  cross  products  of  the  X” 
and  Z"  axis  unit  base  vectors  (direction  cosines). 

cos  ■  cos  Y2  cos  cir^  -  cos  a2  cos  y3 

(26) 

cos  P2  "  cos  a±  cos  Y3  -  cos  cos  03 

(27) 

cos  ^  "  cos  >1  cos  -2  “  009  ®1  cos  Y2 

(28) 

the  unknown  translation  terms  are  determined  from  Equations  17, 
for  point  1  in  Surface  A  (xj[  *  0,  y£  -  0,  -  0)* 

18  and  19 

H”  *  -  cos  CA  -  72  cos  or2  -  *2  cos  ^3 

(29) 

L”  -  -  X2  cos  &2  -  ?!  cos  ?2  -  Z2  005  ^3 

(30) 

KM  «  -  X2  cos  Y^  “  72  cos  Y2  -  z 2  c03  Y3 

(31) 

The  prograi  now  transforms  all  coordinates  in  Surface  A  a^d  Surface  B  to 
the  new  system  by  using  Equations  17.  18  and  19* 
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APPENDIX  D 


CONFUTATION  OF  SURFACE  AREA  OF  INTERNALLY  (ENERATED  SURFACES 


In  Figure  54,  View  J-J  shows  a  view  of  the  surface  of  a  right  elliptical 
cone  between  the  two  arbitrary  cross-sections  indicated  in  the  isometric 
sketch  of  the  cone.  Because  the  program  internal  surface  generator  uses 
the  elliptical  cross -ssction  as  the  basic  generating  element,  elemental 
s**1  face  areas  such  as  A^D  in  Figure  54  are  trapezoids  having,  in 
general,  unequal  nonpar/dlel  sides.  Also,  because  each  elemental  surface 
is  developed  by  equal  elliptical  Parametric  angles,  one  need  compute  the 
area  of  only  cne  elemental  surface  for  each  pair  of  cross-sections 
(providing,  of  course,  the  cross-sections  are  similar). 

The  plane  area  of  trapezoid  ABCD  is 

Aj  »  i  h  (I,  +  L,).  (1) 

In  and  1>2  are  readily  computed  from  the  X,  Y  coordinates  of  points  A,  B, 

C  and  Dj 

h  -V&C  -  xa>2  +  Oc  -  (2) 

h -fit)  -  h'l2  +  <yd  -  *B>2  <3) 

The  trapezoid  height  h  is  computed  indirectly  from  the  projected  area  Ap  of 
the  trapezoid  on  the  AY  plane  in  the  following  manner: 

h  -^hp2+Z2  (4) 

where  h  is  the  projected  length  of  trapezoid  height  h,  and  Z  is  the  dis¬ 
tance  between  cross-sections. 

The  projected  area  of  ADCD  is: 

*p"4hp  <h+I2>-  <»> 

Solving  for  h_  in  Equation  5’ 

P  .  2  AP 

N>  “  (Ij  +  \]  (6) 

The  area  is  computed  from  the  trapezoid  (X,Y)  coordinates  in  tho  fol¬ 
lowing  manner: 
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Noting  triangles  BOD  and  AOC 

Ap  =*  Area  A  BOD  -  Area  AAOC  (7) 

It  is  desirable  for  computational  efficiency  to  compute  A_  using  known  para¬ 
meters.  In  this  case,  sin  6  ^.s  known  and  is  used.  The  p&rametric  equations 
of  the  ellipse  are: 


7  *  b  sin  9  (9) 

where  a  is  the  semi-major  axis,  b  is  the  semi-minor  axis  and  B  is  the  para¬ 
metric  angle.  In  Figure  54,  angle  Sj_  defines  points  A  and  B  and  angle  Bp 
defines  points  C  and  D.  Angle  (J  “  Bp  -  $1* 

The  area  of  triangle  AOC  can  be  computed  by  vector  cross  products. 


From  Equations  8  and  9  for  the  parametric  angles  0i  and  02  defining  points 
A  and  C, 

XA  "  a!  cos  YA  °  sin  (11) 

Xg  "  a^  cos  Bp,  Yq  ■  sin  Bp  (12) 

Substituting  Equations  11  and  12  in  Equation  10, 

AA0C  “  4  (c°&  ®1  sin  -  cos  Bp  sin  0^) 

aA0C  *  4  sin  (02  -  Bj.) 

aAX  “  4  sin  0  (13) 

A  similar  derivation  is  developed  to  obtain  the  area  of  triangle  BOD. 

AB0D  **  4  *2^2  ®  (14) 

Substituting  Equations  13  and  14  in  Equation  7, 

Ap  -  i  sin  8  (a2b2  -  a-^)  (15) 

Substituting  Equation  15  into  Equation  6, 

«  2  [  4  sin  B  (*2b2  -  *1^)] 


Substituting  Equation  16  into  Equation  Uj 


+z2  <i?) 
Finally,  substituting  Equation  17  into  Equation  1, 

u  -  bJih  ,/h 8  (s2b;  -  2  +  ^ 

1  2  VL  h+h  1 

Rearranging  terns, 

Aj  -  j  ^/[s in  B  (ajb2  -  a^)]  2  +  [z  (I,  +  I^)]  2  (18) 

The  total  surface  area  is  computed  by  repeated  evaluation  of  Equation  18 
for  the  particular  surface  generated. 
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Apparon  s 


DERIVATION  C?  CONFIGURATION  FACTOR  TO  A  SPHERE 


The  analytic  jolution  of  the  configuration  factor  to  a  sphere  depends 
upon  the  position  of  t^e  sphere  relative  to  the  plane  of  Stir  face  1.  Inis 
is  clearly  demonstrated  geometrically  by  using  the  Nusseit  unit  sphere 
projection  as  shown  in  Figure  55.  Three  unique  solutions  ore  apparent 
iron  the  five  different,  sphere  positions;  lb  and  Ilia  represent  "limit” 
values  in  each  case.  The  unit  sphere  projections  (crosshatchid  areas) 
shown  in  the  lower  part  of  Figure  55  correspond  respectively  to  sphere 
positions  (not  to  scale)  depicted  in  the  upper  half.  Case  I  results  vher 
the  sphere  is  above  (la)  and/or  touching  (lb)  the  plane  of  SurJUce  1. 

The  locus  on  the  N  wselt  hemisphere  base  is  an  ellipse,  and  varies  froc>  a 
circle  (when  the  spnere  is  vertically  ovei  the  point  0;  to  the  single  tan- 
gency  position  sho’-n  in  To.  When  the  sphere  goes  below  the  plane  of  Sur¬ 
face  1,  the  Case  II  locu3  appears,  and  is  formed  by  the  ellipse  boundary 
on  the  left  and  the  unit  circle  boundary  on  the  right. 

The  locus  projected  on  the  unit  sphere  surface  is  a  circle,  in  every 
case.  The  radius,  b,  of  the  circle  becomes  the  semi-major  axis  of  the 
projected  ellipse.  By  similar  triangles, 

b  „  R 

1  D 

!) 

o  -  “g-  (!) 

where  D  «  distance  from  center  of  sphere  to  origin  of  unit  sphere  and  R  is 
the  sphere  radius.  The  semi-minor  axis,  a,  of  the  ellipse  is  the  projec¬ 
tion  on  the  unit  circle  oi  b;  again  by  similar  triangles. 


JL  -  -ft- 

D  b 


The  area  of  the  ellipse  is 


2  _R2 
D  1? 


•  nab  1 


(2) 


(3) 


The  configuration  factor  for  Case  I  is  the  aroa  of  the  ellipse  divided  by 
the  area  of  the  unit  radius  circle. 
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(0 


The  computation  of  Case  III  proceeds  as  follows,  referring  to  Figure  56, 
a  detailed  drawing  of  the  upper  half  of  the  Case  III  b  locus,  the  area  to 
be  determined  is  Area  ACB. 


By  inspection. 


AACB  "  AAOB  ■  A0'0B  •  %>’* 


Area  AOB  Is  a  sector  of  the  unit  circle,  R  ■  1, 


AOB 


“  h HI  - \  6 


X 


(5) 

(6) 


Given  h,  the  distance  between  centers  0  and  G1,  and  Yq,  the  value  of  y  at 
the  tangency  point  between  the  ellipse  end  the  circle,  the  a*  ?i  of  the 
triangle  O' O'*  1$ 


Votf 


yc 


(7) 


The  area 

<V  =rc) 


CO'B  is  an  elliptical  sector  defined  by  the  parametric  angle  0  at 
the  point  of  tangency. 


CO'B  2 


(B) 


Inserting  Equations  6,  7  and  8  into  Equation  5> 

aacb  “  £  (d  -  hyc  "  (9) 


.'he  unknowns  6,  h,  jc  am  p  must  be  evaluated  in  terms  of  R,  2  and  D.  The 
tangency  point  xr,  y^,  is  determined  as  follows  The  equation  of  fhe  el¬ 
lipse  when  translated  a  distance  h  in  the  x  direction  from  the  origin  0  i3 


(10) 


The  equation  of  the  unit  circle  ia 

x2  +  y2  *  1  (11) 

At  the  intersection  point  the  slopes  are  equal ,  Taking  the  first  derivative 
of  Equation  10, 


Roarranging, 


~y  (x"h)  ^  r  -f?  dy  *  0 

&  =  si  [.a,  oJ  -  j4 

dx  ~y  a-  j  a*-  y 


(12) 


The  slope  at  any  point  x,  y  on  the  circle  i3: 


dv  r,  _  _2S_ 

dx  y 


Equating  Equations  12  and  13, 

i2  (hr*)  «  _  J5- 


Solving  for  h,  at  x 


xc  and  y  "  yc, 

„  (V2;.2) 

h  "xc 


Substituting  Equation  14  irto  Equation  10, 


(13) 


(14) 
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h-*°  liira] 


+  Za--i 

h2 


Reducing  nr/J  rearranging  terms,  and  solving  simultaneously  with  Equation  II 

(“b~)Z  :co2  *  y0z  -  b2  (15) 

2  2 

V  +  yc  -1 

Subtracting  Equation  11  from  Equation  15, 

xeJ  [(f-)2  -  l]  -b2  -  1 


Solving  for  r  , 


(«)2-  i 


X  «  b  /  ft-r-KL 

'  V  b2  -  ,2 


Solving  for  yQ  in  Equation  11, 

y0  -  V 1  -  *c2 

Substituting  Equation  1  and  2  into  Equation  17; 

xc  "J-frtj? 

Substituting  Equation  1,  2  and  16  into  Equation  18, 


./HZ' 

V  ?:2  J7 


The  angle  0  may  be  defined  as 

y 

0  -  tan-1  jB.  -  tan”1  r*° - r- 

K  Vb2  -  yQ2 

Substituting  Equation  1  and  20  into  Equation  21, 


(21) 


The  angl.  0  nay  be  defined  as 


6  •*  tan"  ^ 

r„ 


(23) 


Substituting  Equations  1?  and  20  into  Equation  23, 


y  -  tan 


,  fz2~  1% 

'  Vim; 


n2 

r 


Subrtituting  Equations  1,  2  and  19  into  Bq'u.tion  14, 


D2 


l  -  Z2  [WZ 

1)2  V  Z2  - 


£ 

D2 


(?4) 


(25) 


The  parameter  hy,  in  Equation  9  is  evaluated  by  multiplying  &yj:>tir.n*  20 
and  25, 


hy„ 


hy. 


1 1 


£=z£  /ETif  1 

D2  V  Z2  -  1)2  J 
Vd2  -  r2  •  y/?2  -  z2' 

D2 

w  -  (D2  -  n2s  /  R2  -  z2 
c  D2  V  J52  -  R2 


Z2-  R2 


Let 

and 


9-  «  tr  -  R 

0 


2 


(26) 

'2?j 


(28) 


/E 3T 

▼  O' 

Substituting  Equations  2?  and  28  Into  Equati  *ns  22,  21+  and  26,  respectively, 

0  »  tan-1^-  3  (29) 


0  ”  tan--1  6 


(30) 


and 


hy, 


c  d2 


(31) 


Final!  ,  substitute  Equations  29,  30  aixd  31  into  9!  the  actual  area 
project*;  ox?  the  hemisphere  base  is  tv.dee  Area  ACB  and  ic  divided  h/  the 
base  are  zo  yield  the  configuration  .factor. 
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::  2) 


where  ot  « 

3  *%/p2  "-Z-  »  Z  «  0  and  Z2  * 

By  inspection  of  Figure  55>  tho  projected  *rea  for  -.e  I  is  the  st.m  of  the 

t.'lli Due  evaluated  by  the  Cane  I  formula  anti  the  cri  jf.t  shaped  are;.  r.t  tr- 
mined  by  the  Case  III  ioraula. 


’II 


ei  +  'hi-  0 2  *  * 


In  suauiAry,  referring  to  Figure  55 $ 


<?-  u.  Z  >  r 


D3 


where 


i  -  sl-  ^l4*"3'^)-  /] 


Z  <  0,  Z2  2  R2 


/•  «  /t  j.  /?  n  7  <  p 

tfu  ei  iip  ** 
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